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ABSTRACT 


Using the JVLA, we explored the Galactic center (GC) with a resolution of 0.05" at 33.0 and 44.6 GHz. We 
detected 64 hyper-compact radio sources (HCRs) in the central parsec. The dense group of HCRs can be divided 
into three spectral types: 38 steep-spectrum (ao < —0.5) sources; 10 flat-spectrum (—0.5 < a < 0.2) sources; and 
17 inverted-spectrum sources having o > 0.2, assuming S x v^. The steep-spectrum HCRs are likely represent 
a population of massive stellar remnants associated with nonthermal compact radio sources powered by neutron 
stars and stellar black holes. The surface-density distribution of the HCRs as function of radial distance (R) from 
Sgr A* can be described as a steep power-law X(R) x RT, with T = 1.6 +0.2, along with presence of a localized 
order-of-magnitude enhancement in the range 0.1-0.3 pc. The steeper profile of the HCRs relative to that of the 
central cluster might result from the concentration massive stellar remnants by mass segregation at the GC. The GC 
magnetar SGR J1745-2900 belongs to the inverted-spectrum sub-sample. We find that two spectral components 
present in the averaged radio spectrum of SGR J1745-2900, separated at v ~ 30 GHz. The centimeter-component 
is fitted to a power-law with a, 2 —1.5 30.6. The enhanced millimeter-component shows a rising spectrum 
Qmm = 1.10.2. Based on the ALMA observations at 225 GHz, we find that the GC magnetar is highly variable 
on a day-to-day time scale, showing variations up to a factor of 6. Further JVLA and ALMA observations of the 
variability, spectrum, and polarization of the HCRs are critical for determining whether they are associated with 
stellar remnants. 


Unified Astronomy Thesaurus concepts: Center of the Milky Way; [Galactic center (565)]; Interstellar medium 
(847); Radio continuum emission (1340); Black holes (162); Pulsars (1306); Magnetars (992); Neutron stars 


(1108); Discrete radio sources (389); Radio transient sources (1358); Radio interferometry (1346) 


1. INTRODUCTION 


One of the outstanding open questions that has challenged 
astronomers for many years is the “missing pulsar problem”: 
there are far fewer pulsars found toward the Galactic center 
(GC) than we could expect, given the formation rate of mas- 
sive stars in the central molecular zone of the Galaxy implied 
by the relative abundance of massive stars produced at the GC 


over the past ten million years (e.g., Dong, Wang & Morris. 
2012} [Lu et al. |2013] (Clark et al. |2021). This can in part 


be ascribed to the large foreground scatter-broadening at longer 
radio wavelengths toward the GC, which can lead to a large 
enough pulse broadening that the pulses become indistinguish- 


able (Lazio & Cordes 11998). Several other reasons also com- 


plicate the discovery of GC pulsars, as detailed by 
(2021). However, the discovery of a magnetar associated 


with SGR J1745-2900, located just 3" from Sgr A* 
Rea et al. 2013), indicates that 
the effect of the scattering screen could be up to three orders of 
magnitude smaller than had previously been expected 
‘Bower et al. 2014). Consequently, the question re- 
mains: why haven't more pulsars been seen toward the Galac- 
tic center? Because massive stars clearly form in abundance 
at the GC, and because dynamical friction should cause the 
more massive stellar remnants to be concentrated there, neu- 
tron stars should be abundant and continuously produced in the 


GC region (R ~ 0.5 pc) (Bahcall & Wolf [1976 
Miralda-Escudé & Gould [000] Prah! & Loeb [2004] Alexan- 
der & Hopman 009] Merritt D010j Antonini & Merritt [2012 
| 


One obvious answer to this question of where the pulsars 
are is that the number of “windows” in the scattering screen is 
quite small, so that most pulsars are still too scatter-broadened 


for their pulses to be detected at the wavelengths searched. An- 
other, perhaps more interesting answer is that massive stars that 
form out of the rather highly magnetized interstellar medium 
of the GC tend to themselves be rather strongly 
magnetized, and therefore leave strongly magnetized neutron 
star remnants. That is, pulsars formed near the GC could fre- 
quently be magnetars, which have short lifetimes as recogniz- 
able pulsars (~ 10? — 10? yrs) because of their rapid spin-down 
Beloborodov |2017). Such short lifetimes would limit the num- 
ber of pulsars that could be detected at any one moment to a 
small number, although they could remain detectable as com- 
pact radio sources. Radio continuum surveys of point sources 
can help distinguish these possibilities. We recently published a 
5.5-GHz survey of GC compact radio sources (GCCRs) within 
a radius of ~7.5 arcmin (17 pc) of Sgr A* (Zhao, Morris &| 
[Goss 2020), and concluded that, of the 110 new compact ra- 
dio sources observed down to a 10-c sensitivity limit of 70 uJy, 
most of them fall within the high flux density tail of normal pul- 
sars at the GC (our effort to decrease the 5.5-GHz flux density 
limit with existing, additional data is in progress). Of course, 
there are several other possible assignations for these sources; 
82 of them are variable or transient and 42 have possible X-ray 
counterparts. 

Limited by the VLA angular resolution and confusion from 
the HII continuum emission from Sgr A West, the 5.5-GHz sur- 
vey focussed primarily on regions lying beyond a radius of R ~ 
1 pc from Sgr A*, that is, on regions outside the circumnuclear 
disk. To take the next step in addressing the pulsar puzzle, we 
have recently surveyed the central +0.5 parsecs (+13") around 
Sgr A* at higher frequencies, using existing JVLA Ka and Q- 
band data and X-band observations in the A-array. The high- 
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resolution JVLA observations at 33 and 44.6 GHz were used to 
search for hyper-compact (< 0.1") radio sources (HCRs) in Sgr 
A West, and to study their radio properties and distribution near 
Sgr A*. The motivation for going to higher frequencies in the 
context of constraining the magnetar population comes from the 
discovery by{Torne et al. |(2017) that the spectrum of the mag- 
netar near Sgr A*, SGR J1745-2900, rises at higher frequencies 
to a millimeter/submillimeter plateau. Another magnetar, 1E 
1547.0—5408, has also been seen to display a spectrum rising 
at millimeter wavelengths (Chu et al. 2021). If such a rising 
spectrum happens to be a general characteristic of magnetars, 
then this feature can be used to identify magnetar candidates 
with higher frequency observations, even if radio pulses are not 
detectable. 

This paper is organized as follows: Section 2 describes the 
JVLA observations, data reduction and imaging procedures used 
for identifying HCRs within the central parsec. We also de- 
scribe there our procedure for data reduction and imaging using 
archival ALMA data for measurements of SGR J1745-2900. 
Section 3 presents a catalog of the HCRs found within Sgr A 
West. Three selected cases of HCRs are described in Section 
4, including detailed results on the radio spectrum and variabil- 
ity of the GC magnetar SGR J1745-2900, based on data from 
this paper and from prior publications. Section 5 shows a sta- 
tistical analysis of HCRs by dividing them into three spectral 
types. Possible origins of the HCRs as well as massive stellar 
remnants as a consequence of the mass segregation in the cen- 
tral parsec are also discussed in section 5. Finally, section 6 
summarizes our conclusions. 


2. OBSERVATIONS, DATA & IMAGING 


Following the 5.5 GHz JVLA A-array survey of compact ra- 
dio sources within a radius of 7.5 arcmin (17 pc) in the Galactic 
center region, we have carried out a search for compact radio 
sources within a radius of 13" (0.5 pc) based on the existing 
JVLA high-resolution data at 44.6 GHz (Q-band) and 33.0 GHz 
(Ka-band) as well as our recent JVLA A-array observations at 
9 GHz. The high-resolution and sensitive VLA observations at 
high radio frequencies are crutial in detections of HCRs at a 
level of 100 uJy to a few mJy in the vicinity of Sgr A*. 


2.1. JVLA data, calibration & imaging 


New JVLA observations in the A configuration were carried 
out on 2019-9-21, 2020-11-20 and 2020-12-4 at 9 GHz, in ad- 
dition to the observation on 2014-04-17 that was used to image 
the Sgr A West filament (Morris et al. [2017). The X-band ob- 
servations were all carried out with the VLA standard correla- 
tor setup for wideband continuum covering 2 GHz bandwidth, 
produced from the 8-bit sampler. We also acquired the higher- 
resolution NRAO archival data observed with the JVLA in the 
A-array at 33 and 44.6 GHz on 2015-9-16 and 2015-9-11, re- 
spectively. The Q- and Ka-band observations were using the 
3-bit sampler, producing 64 subbands and covering a total of 8 
GHz bandwidth. All the observations were pointed at a sky po- 
sitioq!] very near Sgr A*. Table 1 summarizes the six sets of uv 
data (columns 1 - 8). The data reduction was carried out using 
the CASAP| software package of the NRAO. The standard cal- 
ibration procedure for JVLA continuum data was applied. 3C 
286 (71331-43030) was used for corrections for delay, bandpass 
and flux-density scale. 


! RA(2000) = 17:45:40.0383, Dec(J2000) = —29:00:28.069 
? http://casa.nrao.edu 


J1733-1304 (NRAO 530) and J1744-3116 were used for com- 
plex gain calibrations. In addition, corrections for the time vari- 
ation of the bandpass across each baseband due to residual de- 
lays were determined using NRAO 530 based on the model dis- 
cussed in [Zhao, Morris & Goss |(2019) . The accuracy of the 
flux-density scale at the JVLA is 396—596, limited by the uncer- 
tainty of the flux density of the primary calibrator, Cygnus A 


(Perley & Butler 12017). Following the procedure for high dy- 
namic range (DR) imaging that we developed recently 


and applying it to the Sgr A data with 
CASA, we further corrected for the residual errors in phase. 
The integration time of the calibrated Q and Ka-band data was 
averaged into 30-sec time bins, so that the intensity loss at the 
13" outer radius of the search region is less than 1096 due to 
the smearing effect caused by the Earth's spin. After correcting 
for the residual delay, we also binned the spectral channel data 
to channel widths of 32 and 16 MHz for the Q- and Ka-band 
data, respectively, to ensure that the intensity loss for a point 
source caused by bandwidth smearing is less than 10% at the 
radius of 13". The short baseline data (< 500kA) were filtered 
out in order to avoid contamination by the extended emission 
in Sgr A West on scales > 0.4". The Ka-band image shown in 
Fig. 1 is made with CASA task CLEAN from the baseline data 
between 500-4,000 kA, achieving an rms of 8 uJy beam ! with 
a FWHM beam of 0.079" x0.031 (—11°). This high-resolution 
image shows numerous compact radio sources in the central 
parsec region. 

For the X-band data, after corrections for the residual errors, 
the visibility data were averaged to a time bin of 16 sec while 
the original channel width of 2 MHz was kept. We made im- 
ages at 9 GHz using the multi-frequency synthesis algorithm 
(MFS: [Rau & Cornwell [2011) with the 1024 spectral channels 
covering the 2 GHz bandwidth. Also we filtered out the short 
baselines (« 100 kA), and constructed two 9-GHz images for 
both the 2019 and 2020 epochs. The second epoch image was 
made with the two observations on 2020-11-20 and 2020-12- 
04. Hereafter, we use the mean epoch, 2020-11-27, for this 
image. The rms noise for the 2020-11-27 X-band image is 5 
uJy beam !. We reconstructed the 2014-04-17 image with the 
calibrated X-band A-array data (Morris et al. |2017) by using 
only the long-baseline data (> 100 kA). The rms noise for the 
resulting 2014-04-17 X-band image is 4 pJy beam !. 

The specified parameters for the final images at Q-, Ka- and 
X-bands are summarized in Table 1 from columns 9 to 11. 


2.2. ALMA data, calibration & imaging 


We acquired archival data from the Atacama Large Millime- 
ter Array (ALMA), observed by [Tsuboi et al. (2019) at 225.75 
GHz in Cycle 5 (2017.1.00503.S). Following the ALMA CASA 
guide for Cycle 5 data reduction, we executed the pipeline script 
"scriptForPI.py" to produce calibrated ALMA data in CASA 
Measurement Set format. The ALMA datasets are composed 
of ten 1-hr observations of IRS 13E in the array configuration 
of C43-10 within a two-week interval between 2017-10-6 and 
2017-10-20. The observing field covers Sgr A* and the mag- 
netar SGR J1745-2900 in addition to IRS 13E. The pipelined 
ALMA images appear to be marred by severe residual errors. 
Both IRS 13E and SGRA J1745-2900 were buried in the side- 
lobes and artifacts produced by the residual dirty beam. Then, 
following the recipe for dynamic range imaging with wideband 


data (Zhao, Morris & Goss |2019), we corrected the residual 
errors by utilizing the antenna-based closure relations (Thomp-| 
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TABLE 1 


LoG oF JVLA DATASETS AND IMAGES 


UV data Images 
Project ID Array Band v Av At HA range Epoch Weight (@maj-Omin, PA) rms 
(GHz) (GHz) (sec) (day) (R) (arcsec, arcsec, deg) (Jy beam!) 

(1) (2) (3) (4) (5) © (7) (8) (9) (10) (11) 
15A-293 A Qt 44.6 8 3 0.6 — +2!.5 2015-09-16 0 0.078, 0.032, 12 17 

Ka! 33.0 8 2 -2h8— «0.3 2015-09-11 —0.3 0.079, 0.031, —11 8 
14A-346 A xi 9.0 2 2. -3h4— «3^3 2014-04-17 0 0.36, 0.15, —6 4 
19B-289 A xi 9.0 2 2. 40^.6— 42^.6 2019-09-21 0 0.36, 0.15, —6 7 
20B-203 A xt 9.0 2 2. -0'.3 —+2'.9 2020-11-20 0 0.36, 0.15, -6 5 

-1^.5 — 414.7 2020-12-04 


(1) JVLA program code of PI: Mark Morris for 19B-289 and 20B-203; PI: Farhad Yusef-Zadeh for 14A-346 and 15A-293. (2) Array configurations. (3) JVLA band code; "X","Ka", and "Q" stand for the VLA bands in the 
ranges of 8.0 — 12.0 GHz, 26.5 — 40.0 GHz and 40.0 — 50.0 GHz (https://science.nrao.edu/facilities/vla/docs/manuals/oss2013B/performance/bands). (4) observing frequencies at the observing band center. (5) bandwidth. (6) 
integration time. (7) hour angle (HA) range for the data. (8) date corresponding to the image epoch. (9) robustness weight parameter. (10) FWHM of the synthesized beam. (11) rms noise of the image. t Correlator setup: 64 


channels in each of 64 subbands with channel width of 2 MHz f Correlator setup: 64 channels in each of 16 subbands with channel width of 2 MHz. 


4 Zhao, Morris, Goss 


son, Moran & Swenson |2017) and constructed a refined ALMA 


image from the data observed at the first epoch on 2017-10-6 
as a trial case using CASA task tCLEAN with the MFS algo- 
rithm and robust weighting with Briggs parameter R=0 
[|[1995). An rms noise of ø = 20 uJy beam ! was achieved, 
with a FWHM beam of 0.024" x0.017" (82°). The magne- 
tar SGR J1745-2900 was successfully detected with a S/N ra- 
tio of 70. We then made 2D-Gaussian fits to both Sgr A* 
and the magnetar, finding flux densities of 3.0312-0.013 Jy and 
1.41+0.08 mJy at 225 GHz for Sgr A* and the magnetar, re- 
spectively. Fig. 2a shows the ALMA continuum image of the 
magnetar at 225.75 GHz. Based on the procedure and input pa- 
rameters for the CASA tasks that were used for the trial case, 
we coded the detailed CASA reduction steps for the ALMA 
data into a CASA-Python script. Using this script for correc- 
tions of residual errors, we further processed the datasets from 
all ten epochs observed on 2017-10-6, 2017-10-7, 2017-10-9, 
2017-10-10, 2017-10-11, 2017-10-12, 2017-10-14, 2017-10- 
17, 2017-10-18, and 2017-10-20. The flux densities for the 
magnetar were determined by fitting a 2D-Gaussian and are tab- 
ulated in Table 3 (see Section 4). 
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FIG. 2.— (a) Contour image of the magnetar SGR J1745-2900 ob- 
served on 2017-10-6 using ALMA at 225.75 GHz with a FWHM beam of 
0.024" x 0.017" (PA 82°). Contours are Cim x (—5.7, 5.7, 8, 11.3, 16, 22.6, 32, 
45.3 ), where Cim = 22uJy beam"!. (b) ALMA image of the magnetar SGR 
J1745-2900 (top-right: a point source) and IRS 21 (bottom-left: a complex) 
observed on 2019-10-14 at 318 GHz with a FWHM beam of 0.087" x 0.059" 
(89.59). Contours are Cim x (—6, 6, 8, 11, 15), where Cim = 55uJy beam"! . The 
FWHM beams are illustrated at the bottom-left corner in each image. Numbers 
on the grey scale bar are in units of Jy beam"! . 


We also re-processed ALMA archival data (2015.1.01080.S) 
for observations at 343.49 GHz by {Tsuboi et al. |(2017). The 
observations were carried out at four epochs 2016-4-23, 2016- 
8-30, 2016-8-31 and 2016-9-8, for durations of 3h, 1h, 2h and 
3h, respectively. The 2016-4-23 observation was in the C36- 
2/3 array configuration, and other three observations were in 
the C40-6 configuration. The pipeline-calibrated datasets were 


obtained by executing the pipeline script "scriptForPI.py". We 
then adjusted the input parameters in the CASA-Python script 
used for imaging the 225-GHz data, and applied the script to 
the 343-GHz data. The spatial resolution of the first 343-GHz 
image from epoch 2016-4-23 is relatively poor, with FWHM 
beam of 0.35" x 0.33" (—79?), and the emission from the magne- 
tar SGR J1745-2900 appears to be contaminated by surround- 
ing extended emission. The final image of SGR J1745-2900 
was made by applying a high-pass baseline filter, so that only 
long-baseline data were included (> 100 kA), thus filtering out 
extended ( 2") emission features. An rms noise of 13 mJy 
beam ! was achieved for the epoch 2016-4-23, and the flux den- 
sity of 3.29+0.34 mJy was determined for the magnetar at 343 
GHz. The observations from the later epochs at 343 GHz have 
a typical angular resolution of 0.1". The rms noise values were 
0.15, 0.05 and 0.12 mJy beam ! for the 2017-8-30, 2017-8-31 
and 2017-9-8 images, respectively. The magnetar SGR J1745- 
2900 is detected in all four epochs at 343-GHz, and the flux 
densities are determined at a level of 10c or better using 2D- 
Gaussian fitting. The measurements are reported in Table 3, 
including a 596 error in the flux-density calibration 
al. 2018), 

A high-resolution (0.087" x0.059", PA = 89.5?) ALMA ob- 
servation was carried out at 320 GHz in the C43-6 array con- 
figuration on 2019-10-14 (2018.A.00052.S of P.I. Mark Mor- 
ris). The data was initially processed with the pipeline script 
"scriptForPI.py". Subsequently, we made further corrections 
for the residual errors with the CASA-Python script described 
above and then imaged the region containing the GC magnetar, 
Sgr A*, and IRS 21 with the two lower-frequency sub-bands 
centered at 318 GHz, which have relatively stable phase and 
less contaminations from molecular lines in the circumnuclear 
disk (CND). An rms noise of 55 pJy beam ! was achieved. 
The magnetar is significantly detected with a flux density of 
1.32+0.15 mJ y. The uncertainties include a 5% error in the de- 
termination of the flux-density scale for ALMA observations. 
Fig. 2b shows the 318-GHz image of the region containing 
SGR 71745-2900 and IRS 21. 


3. HYPERCOMPACT RADIO SOURCES 
3.1. Search procedure and selection criteria 


The angular resolutions of the A-array observations at Ka- 
and Q-bands, 0.05", are nearly two orders of magnitude smaller 
in beam area than the 5.5-GHz beam used in our previous study 
of the GCCRs (Zhao, Morris & Goss |2020). At such a high 
spatial resolution, most of the emission components in the Sgr 
A West HII region have been resolved out. The emission from 
the HII sources produced overwhelming confusion at 5.5 GHz, 
which was the main issue preventing us from unambiguously 
detecting GCCRs in Sgr A West, given the relatively low angu- 
lar resolution of 0.5". The Ka-band JVLA, with both its supe- 
rior angular resolution and sensitivity (~ 8 Jy beam ), is well 
suited for significant detection of point-like hyper-compact ra- 
dio sources (HCRs) at a level of >0.1 mJy. Thus, we have 
an unprecedented opportunity for study of the compact radio 
sources associated with massive stellar remnants, such as pul- 
sars, magnetars and accreting compact stellar remnants. Unlike 
the free-free emission of the HII components, the radio emis- 
sion from such objects is expected to be nonthermal. Our search 
has therefore been focused on using the Ka- and Q-band im- 
ages to identify nonthermal HCRs within the Galaxy’s central 
parsec. We proceeded in three steps as follows: 
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Fic. 3.— Images of individual HCRs, color for 44.6 GHz intensity and contours for 33 GHz from an image made with a long baseline filter of > 500 kA except 
for HCR39 which was made with a longer baseline filter > 1000 kA to separate the point-like source from the contamination of the extended emission associated 
with the radio core IRS 13E. The IR stars (E2, E3, and E4) detected in the IRS 13E region are marked with white dots (see HCR38 and HCR39 
panels). The color wedges scale the intensity of 44.6-GHz emission in units of Jy beam™!. The contours are oim x 2"/? where n starts from n; until reaching the 
source peaks. The values of cim and n, corresponding to the individual HCR image panels are listed in Column 11 of Table 2. Bottom-right corners are the FWHM 
beams at 44.6 and 33 GHz: 0.08" x 0.03" (PA = 12°) in blue and 0.08" x 0.03" (PA = -119) in white, respectively. 


(1) We initially found approximately a thousand compact ra- 
dio sources having angular size of < 0.5" based on a Ka-band 
image constructed including all baselines with a robust weight- 
ing parameter, R=0.25 (Briggs |1995), and a FWHM beam of 
0.12" x0.06" (PA 2 -11?). A sensitivity threshold of S33GHz/Orms 
> 6, Oms © 7 py beam™!, was applied in the initial search. 
Bright intensity maxima in the HII emission components might 
have been included in the compact-source sample, producing 
false detections for the nonthermal compact radio sources. To 
ensure that the ultimate sample contains only compact nonther- 
mal radio sources, two further steps were carried out: 


(2) The Ka-band image was reconstructed with a robust weight 
R=—0.3 to down-weight the contribution from short baselines 


and also a high-pass baseline filter (> 500 kA) was applied. 
So we can separate point-like sources from extended emission. 
The final cleaned image was convolved with a beam having 
FWHM = 0.08" x0.03" (PA = —11?), similar to the synthesized 
beam of the Q-band observations, but keeping the position an- 
gle of the Ka-band synthesized beam determined by the uv- 
data sampling. We narrowed the list for those sources that 
are detected at Ka-band by keeping only those sources above 
a S33/0rms = 15 threshold and limiting the size of the major axis 
to be less than 0.1 arcsec (Omaj < 0.1") based on a 2D Gaus- 
sian fitting. The number of candidates in the list was thereby 
reduced to ~100. The remaining sub-sample only contains the 
hyper-compact members of the initial sample. 


6 Zhao, Morris, Goss 


(3) Finally, we used the high-resolution Q-band image, with 
Cms = 17 Jy beam !, FWHM = 0.08" x0.03" (PA = 12°), to 
narrow down the list of HCRs by imposing a conservative limit 
of S44.6/cim > 10 on the detection significance and an upper 
limit on the positional offset between Ka- and Q-band positions 
less than 3c, where c is given in Table 2 with a typical value of 
a few milli-arcsec (mas), based on 2D-Gaussian fits at the two 
frequencies. We note that the locations of HII peaks depend 
on the uv-sampling. The uv data in the Q-band observations 
were primarily sampled in a positive hour-angle range while 
the Ka-band data were sampled in a negative hour-angle range. 
Consequently, the sources having a significant offset between 
the Q- and Ka-band images are suspected to be HII peaks and 
are therefore rejected from the HCR sample. 

The equal beam sizes of the Q- and Ka-band images allow us 
to determine reliable spectral indices that will be used to further 
distinguish the source types. We do not use spectral indices as 
a selection parameter for the HCRs, given a wide range of val- 
ues for the spectral indices covered by the radio sources associ- 
ated with pulsars, magnetars and stellar-mass black holes. We 
note that the primary 15c significance cutoff used for Ka-band 
sources and the lower 10c cutoff for Q-band were chosen for 
the derivation of significant spectral indices, given the poorer 
sensitivity of the Q-band data. Neglecting fitting errors, for a 
weak point-like source, we expect a maximum uncertainty of 
Oa ~ 0.4 in spectral index. 

Those candidates having consistent results from the Gaussian 
fiting to the Ka- and Q-band images were included in the final 
sample of HCRs, consisting of 64 members. We note that the 
conservative search criteria for HCRs may miss highly variable 
and relatively weaker sources having a steep spectrum, given 
that the sensitivity of Q-band observations is a factor of 2 poorer 
than that of Ka-band data. For example, a hyper-compact radio 
source in the IRS 7SW region (see Fig. 1) was rejected from 
the HCR sample because its Ka-band flux density below the 
15c threshold. 


In short, from the JVLA A-array images observed at 33 and 
44.6 GHz on 2021-9-11 and 2021-9-11, we have identified 64 
hyper-compact radio sources (HCRs) located inside Sgr A West 
within a radius of 13" from Sgr A* based on their compactness 
(a size of Omaj < 0.1") at a conservative significance level of 
S/o > 15 at 33 GHz and S/o > 10 at 44.6 GHz. Three excep- 
tions which were not detected at 44.6 GHz on 2015-9-16 are 
also included: HCR22, HCR32 and HCR64. HCR22 appears to 
be a strong candidate for the object that powers the X-ray PWN 
G359.95-0.04 (see Sec.4.1). HCR64 is a compact transient ra- 
dio source associated with the microquasar discovered during 
the 2005 outburst of the X-ray transient XJ174540.0290031 
[Muno et al. 
Zhao et al. {2009}; further discussion of this source is given in 
Sec.4.2. HCR32 appears to be a highly variable source. It was 
detected with a flux density of 0.50.1 mJy at 9 GHz at the 
epoch of 2014-4-17 but no significant detections were made in 
the 2019 and 2020 epochs at 9 GHz. 

One of the sources, HCR49, is associated with the high- 
velocity head-tail radio/IR source designated the "Bullet" (Yusef- 
Zhao etal. [2009]: it 


now appears to consist of at least two compact components. 


3.2. Images of hyper-compact radio sources 


Fig. 3 shows the images of all 64 HCRs in individual panels 
labelled with the HCR identification numbers, the colored im- 


ages represent the 44.6-GHz HCRs and the superimposed con- 
tours are for the 33-GHz images made with baselines exceeding 
500kA. For those panels containing more than one HCRs, each 
of the HCRs is labelled with its HCR ID. The green plus sym- 
bols mark the positions determined from the 44.6-GHz data. 
For HCR22, HCR32 and HCR64, which were not detected at 
44.6 GHz, the 33-GHz positions are shown. 


3.3. A catalog for hyper-compact radio sources 


Table 2 lists 64 HCRs. Column 1 gives the ID numbers 
for the members of the HCR catalog. Columns 2 and 3 are 
the J2000 coordinates, omitting the common part of 1745” in 
Right Ascension (RA) and —29?00' in Declination (Dec). With 
the exception of HCR22, HCR32, and HCR64, the Q-band co- 
ordinates of the HCRs were determined with respect to Sgr A*'s 
equitorial coordinates at the J2000 epoch 
2014). The positions of non-Q-band sources are determined 
using the Ka-band data. The lo uncertainties in the last dig- 
its of RA and Dec. are given in parentheses. We note that, 
throughout this paper, the positions of the HCRs in the figures 
are labelled in the J2000 equatorial coordinate system with re- 
spect to the position of Sgr A*. Columns 4 and 5 list the flux 
densities at 44.6 and 33 GHz with 1c uncertainties. Columns 6, 
7 and 8 show the source sizes, deconvolved from the telescope 
beam, including major (Omaj) and minor (Onin) axes, as well as 
the position angle (PA), all with their 1c uncertainties. Column 
9 lists the spectral index, o (S ~ v^), determined from the flux 
densities at 44.6 and 33 GHz. The parameters used to plot the 
contours in Fig. 3 for individual HCRs are listed in Column 
10: Cim representing the rms variations of the background re- 
gions near the sources and n, being the integer corresponding 
to the multiplicative factor, 29", specifying the lowest contour. 
Column 11 gives a brief note on individual HCRs. 

Finally, we examined the HCRs that can be spatially asso- 
ciated with the 7mm-IR(3.8m) sources 
by transfering their coordinates to the J2000 equatorial 
coordinate system used in this paper, in which Sgr A* is located 
at RA(J2000) = 17:45:40.0409, Dec(J2000) = —29:00:28.118 
(Reid & Brunthaler |2014). After corrections for the system er- 
rors (ARA = 0.0026*, ADec = —0.049") that mainly caused by 
the difference of the position of Sgr A* between|Reid & Brun-| 
and (2014), 22 HCRS are 
identified with possible 7mm-IR sources within an uncertainty 
of 25 mas, the typical positional accuracy given for the 7mm-IR 
sources. The ID numbers of the possible 7mm-IR counterparts 
are given in Column 11 as well. 


4. TYPICAL CASES FOR HYPER-COMPACT RADIO SOURCES 
4.1. HCR22 and the X-ray PWN G359.95-0.04 


HCR22 is located at the northern termination of the extended 


X-ray source G359.95-0.04 
suggested to be a pulsar wind nebula (PWN) (Wang et al. [2006] 
[Muno et al. |2008). The source at 33 GHz can be character- 
ized by a head (HCR22) with a tail extending ~ 0.3" toward 
the south (Fig. 4a). The radio morphology is consistent with 
the X-ray structure, including the immediate angle of the tail, 
but with a much smaller scale. The radio emssion was also de- 
tected at 9 GHz in the three epoch images 2014-4-17 
lal. 2017), 2019-9-21 and 2020-11-27 (this paper) with a larger 
FWHM of 0.36" x0.15" (-6°). At 9 GHz, the source shows no 
signficant variation in flux density, with measured flux densi- 
ties of 1.70+0.09 mJy, 1.53+0.06 mJy and 1.693-0.05 mJy at 


(a) 33 GHz (b) 9 GHz 
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FIG. 4.— (a) The Ka-band image made with the JVLA A-array data using all 
available baselines, showing the radio structure of HCR22, a candidate source 
of particles and energy powering the PWN. The contours are 5 uJy beam"! x 
(4, 4, 5.7, 8, 11.3, 16) with FWMH = 0.11"x0.04", —11? (bottom-right). 
(b) JVLA 9-GHz images at epochs 2014-4-17 (color) and 2020-11-27 (con- 
tours). The contours are 50 uJy beam ! x (64, 80, 96, 112, 128) with FWMH 
= 0.36" x 0.15", —6? (bottom-right). The green plus marks the peak position at 
33 GHz. The three black ellipses mark the peak positions of the counterpart at 
9 GHz at the three observing epochs as labelled. The size of each ellipse scales 
the 1 c uncertainty of the peak position determined from a 2D-Gaussian fit- 
ting. The color wedge scales the intensity in units of Jy beam! for the JVLA 
images. (c) The Chandra X-ray image (angular resolution = 0.5", from[Zhu et] 


2018) of the PWN G359.95-0.04 2006). The green box at the 


tip of the PWN marks the region of HCR22 as imaged at 9 GHz in panel b. 


the three epochs, respectively. A spectral index of —1.0 + 0.1 
is determined by the flux density of 0.43--0.03 mJy, integrated 
over the entire source, at 33 GHz and the mean flux density of 
1.64+0.06 mJy averaged by the three epochs’ data at 9 GHz. 
Assuming 033/44.6 ^» —1, a peak intensity of < 60 wy beam! 
at 44.6 GHz is extrapolated from the 33 GHz image (see Fig. 
4). The non-detection of HCR22 at Q-band is consistent with 
its steep spectrum. In addition, we noticed that the peak posi- 
tion of the source at 9 GHz moved towards north as time in- 
creases. From a least-squares fitting of the source position at 
the three epochs, we find a significant proper motion of the 
compact source at 9 GHz in the declination direction: giving 
Ha 2 —2.0 + 1.0 mas y7! and us = 7.3 + 1.0 mas y! (Fig. 4a). 
For a distance of 8 kpc to the Galactic center, this proper motion 
corresponds to a velocity of 270 +40 km s^! projected onto the 
plane of the sky, which is consistent with the source's projected 
proximity to Sgr A*. The location, steep spectrum, head-tail 
structure and orientation of HCR22 at 33 GHz along with the 
significant northward motion of the 9-GHz peak, suggest that 
HCR22 is plausibly the candidate source of energetic particles 
that are responsible for the X-ray emission of PWN G359.95- 


0.04. 


4.2. Microquasar of X-ray transient CXOGC 
J174540.0-290031 


HCR64 is the radio counterpart of the microquasar associated 
with the X-ray transient CXOGC J174540.0-290031 that was 


discovered by the Chandra X-ray observatory 
2005 2005). The radio emission from this tran- 


sient was found with the VLA 
[2009). Fig. 5 shows the 2015 Ka-band image (contours) over- 
laid on the color K-band image of 2005. A double source was 
detected in the 2005 VLA observations at 22.5 GHz 
[[2009) with a SW compact component (K40) associated with 
the core at the flux density of 3.443-0.1 mJy. The position offset 
between HCR64 (2015) and the SW component, K40, observed 
in 2005 appears to be insignificant, so we follow 
in identifying this as the core. The bright component 
0.5" NE of the core in the 2005 image is not detected in 2015. 
The NE component was 3.5-£0.1 mJy during the 2005 outburst, 


which appeared to be launched from the microquasar (Muno| 


2005, 2005 2005 
2009). The core (K40) was detected with a flux density 


of 1.12-0.1 mJy at 22.5 GHz at the early epochs in 1991 and 
1999 (Zhao et al. |2009), representing quiescent level of the 
microquasar. 

In our Ka-band image, which was observed 10 years after the 
2005 outburst, the flux density of this compact radio source at 
33 GHz was 0.60+0.07 mJy. The 33-GHz flux density appears 
to be consistent with that observed at 22.5 GHz during 1991 
and 1995, when the microquasar was apparently in its quiescent 
state. From the ratio of the flux densities at 33 and 22.5 GHz, 
we determine a spectral index of a = —1.6+0.4 assuming that 
the flux density in the quiescent state is not variable. The 3- 
c upper limit of 0.15 mJy at 44.6 GHz is consistent with this 
derived steep spectrum. 
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FIG. 5.— Images of the radio emission from the micro-quasar correspoding 
to the X-ray transient CXOGC J174540.0-290031, observed with the JVLA 
in the A-array at 33 GHz (contours) and 22.5 GHz (colors) on 2015-9-11 and 
2005-9-16, respectively. At the bottom-right corner, the FWHM beams of 
0.11" x0.04" (2119) and 0.2" x 0.1" (0?) at 33 and 22.5 GHz are shown with a 
white and filled-blue ellipses. Contour values are Cim x (5.7, 8, 11.3, 16, 22.6 
and 32), where cim = 7uJy beam"!. The two plus symbols mark the positions 
of K40 at 22.5 GHz in 2005 and of HCR22 at 33 GHz in 2015. 
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TABLE 2 
A CATALOG OF HYPERCOMPACT RADIO SOURCES AT 44.6 AND 33 GHZ 


ID RA(J2000) Dec(J2000 [Sc o]46 [Scola  OmajEO Omin tao  PAto @44.6/33 EO Cim, N1 Notes 
(>) C") (mJy) (mJy) (mas) (mas) (deg) (Jy bm) 

a) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 

HCRO1 40.6981(3) 18.323(6) 0.95+0.11 1.00-E0.16 74.419. 56.417. 36.423. -0.17+0.65 7.0, ID9? 

HCRO2 40.5194(4) 27.445(7) 0.20+0.02 0.15-E0.02 68.411. 26.430. 30.445. 0.96+0.60 TQ: 

HCRO3 40.3836(1) 23.150(3) 0.27+0.02 0.20+0.02 57.414. 40.48. 150.436. 1.00 0.48 7.0, 

HCR04 40.2641(1) 33.686(3) 0.26+0.02 0.22+0.01 45.413. 16.48. 100.413. 0.55+0.38 7.0, IRS 9W?! 

HCROS 40.2679(1) 29.243(3) 0.22+0.02 0.19+0.02 unresolved 0.49 + 0.52 7.0, ID4?5 

HCRO06 40.2629(1) 27.217(1) 0.80+0.04 0.75-E0.04 20.46. 15.49. 80.442. 0.2140.25 7.0, ID2?5 

HCRO07 40.1708(1) 29.788(1) 2.54+0.13 1.90+0.10 20.8. 8.48. 56.447. 0.96+0.24 7.0, Magnetar! 

HCRO8 40.1449(4) 16.419(13) 0.23+0.02 0.2140.02 19.410. 13.420. 80.424. 0.30+0.32 7.0, 

HCRO9 40.1286(1) 29.060(4) 0.2140.02 0.2740.02 49.423. 18.411. 84.430. —0.83 0.46 7.0, 1D37> 

HCRI10 40.1166(1) 27.526(1) 0.36+0.03 0.3340.03 15.416. 4.417. 22.444. 0.294+0.36 7.0, IDI? 

HCR11 40.1067(6) 22.277(9) 0.1140.02 0.1740.01 61.431. 26.49. 82.421. -1.45+0.68 7.0, IRS 7E? 

HCR12 39.9905(6) 30.154(12) 0.19+0.01 0.3340.03 51.410. 44.46. 38.418. -1.83+40.42 7.0, 

HCR13 39.9683(6) 29.907(12) 0.25+0.03 0.40+0.04 95.418. 34.46. 94.428. -1.56 0.57 8.5, 

HCR14 39.9502(1) 30.636(3) 0.3140.02 0.40+0.04 74.410. 59.49. 158.464. —0.85 0.36 7:0, 


HCRI5 39.9380(3) 31.162(6) 0.16+0.02 0.29+0.02 72.+16. 18.433. 90.420. -1.97+0.53 8.0, 
HCRI16 39.9350(3) 31.370(6) 0.59+0.07 0.82+0.05 81.417. 62.418. 134.432. —1.09 + 0.43 8.0, 
HCR17 39.9151(3) 28.563(7) 0.52+0.07 0.69+0.11 81.429. 60.421. 166.490 -0.94+ 0.66 7.0, 
HCR18 39.9096(1) 30.895(4) 0.1740.02 0.3140.03 48.414. 29.420. 134.424. —1.99+ 0.50 8.0, 
HCR19 39.9013(1) 30.053(3) 0.41 +0.03 0.54+0.04 82.49. 60.411. 98.424. —0.91 0.34 8.0, 
HCR20 39.8952(6) 28.230(12) 0.20+0.03 0.38+0.04 60.423. 56.423. 144.459. -2.13 +0.61 7.0, 
HCR21 39.8883(3) 31.088(7) 0.32+0.04 0.56+0.03 87.420. 64.416. 158.477. —1.86 + 0.49 8.0, 


Mini-cavity’ 
Mini-cavity” 
e 
Mini-cavity’ 
Mini-cavity?, ID187 


e 


HCR22 39.8408(3) 19.316(6) | «0.06 0.15+0.02 73.416. 50.461. 73.482. N/A 7.0, PWN? 

HCR23 39.8407(1) 29.442(4) 0.92+0.08 1.00+0.11 84.412. 44.430. 14.48. -0.2840.46 8.0, Mini-cavity'?, ID207> 
HCR24 39.8345(1) 29.548(3) 0.19+0.01 0.20+0.01 unresolved -0.1740.34 8.0, Mini-cavity!®, ID19?5 
HCR25 39.8272(1) 29.810(4) 0.88+0.08 1.2940.13 49.418. 34.419. 20.452. -1.2740.45 890, IRS 13E-SE®, ID32>> 
HCR26 39.8247(1) 29.521(3) 0.29+0.02 0.25 0.02 unresolved 0.49+0.42 8.0, IRS 13E-NE!!, 1D21?° 


IRS 13E-SE° 
Mini-cavity 

IRS 13E-NE!!, 1D237° 

IRS 13E-NE!!, 1D227° 
IRS 13E-NEP 
Mini-cavity!* 

IRS 13E-SES, ID307 
Mini-cavity! 

IRS 13E-NE!!, ID245 

IRS 13E-NP, 128? 


HCR27 39.8236(3) 30.039(7) 0.2140.03 0.37+0.04 45.422. 25.428. 133.452. —1.88 + 0.64 8.0, 
HCR28 39.8235(3) 31.825(6) 0.35+0.04 0.45+0.03 73.416. 36.48. 20.416. -0.83+0.47 7.0, 
HCR29 39.8210(1) 29.311(4) 0.30+0.03 0.50+0.03 44.412. 34.420. 148.4336. —1.70 +0.41 8.0, 
HCR30 39.8181(1) 29.475(4) 0.20+0.02 0.28+0.02 58.45. 22.42. 164.43. -1.12+0.47 8.0, 
HCR31 39.8180(1) 28.902(6) 0.63+0.03 0.93+0.06 87.45. 42.47. 98.47. -1.291 026 8.0, 
HCR32 39.8165(2) 30.861( 7) «0.06 0.42+0.03 40.45. 11.422. 72.44. N/A 7.0, 
HCR33 39.8121(1) 29.860(3) 0.35+0.03 0.58+0.05 48.414. 16.410. 14.422. -1.68+0.43 8.0, 
HCR34 39.8107(1) 31.739(4) 0.30+0.03 0.45+0.04 93.414. 45.46. 8.412. -1.35 0.46 7.0, 
HCR35 39.8075(1) 29.300(3) 0.42+0.03 0.434+0.05 55.410. 22.426. 62.424. -0.08+0.42 7.0, 
HCR36 39.8007(3) 29.020(6) 0.31+0.04 0.58+0.04 42.48. 27.416. 56.427. -2.08 + 0.50 8.0, 


-dO 000 -1-):2-1)-]100 1-1 -]1-1]-1OQO00 00000 ON -1-1UA tA — QN NO CA OO MAA 


HCR37 39.7967(6) 30.914(9) 0.18+0.04 0.4140.03 72.425. 41.455. 36.444. -2.7340.79 8.0, IRS 2-N!° 
HCR38 39.7961(1) 29.662(3) 5.50+0.38 3.354042 78.418. 42.418. 52.421. 1.65+40.47 8.0, 12 IRS 13E?, ID33/? 
HCR39 39.7941(3) 29.844(9) 1.2740.12 1.1440.09 76.417. 29.45. 90.47. 0.364041 80,6 IRS 13E-S!7, ID34? 
HCR40 39.7895(1) 28.540(4) 0.58+0.05 0.6940.04 94.49. 54.43. 91.45. -0.58+0.34 8.0, 6 Triplet? 
HCR41 39.7891(3) 29.508(7) 0.90+0.12 1.10+0.07 61.+20. 41.+30. 150.+38. -0.6740.49 8.0,7 IRS 13E-WP? 
HCR42 39.7857(3) 31.239(6) 0.67+40.05 1.00+0.06 84.411. 52.47. 166.412. -1.3340.32 8.0, 7 IRS 2 
HCR43 39.7855(1) 28.421(4) 0.3040.03 0.48+0.05 68.416. 46.411. 161.430. -1.5640.49 80,6 Triplet!® 
HCR44 39.7799(1) 29.587(4) 1.23+0.13 1.234010 46.418. 22.417. 36.433. 0.0040.45 80,7 IRS I3E-WP, ID35? 
HCR45 39.7764(3) 31.938(3) 1.48+0.09 1.15+0.08 96.410. 40.44. L+5. 0.8440.31 8.0, 8 IRS 2L, ID36” 
HCR46 39.7721(1) 29.678(7) 0.64+0.06 0.68+0.04 90.411. 56.414. 78.416. -0.20+40.36 8.0, 7 IRS 13E-WP? 
HCR47 39.7721(2) 31.154(9) 0.18+0.03 0.15+0.03 44.418. 21.6. $82.10. 0.61+40.90 7.0, 5 IRS 2 
HCR48 39.7699(1)  29.479(7) 0.10+0.02 0.19+0.01 35.424. 16.414. 59.434. -2.1340.73 8.0, 7 IRS 13E-W!? 
HCR49 39.7697(3) 25.913(7) 0.17+0.02 0.2140.04 65.419. 39.413. 78.414. -0.7040.78 7.0, 7 Bullet SE* 
E 39.7671(4) 25.909(6) 0.13+0.02 0.16+0.04 63.414. 23.405. 5.10. -0.69+41.00 M Bullet NW* 
HCR50 39.7677(1) 29.824(3) 0.27+40.02 0.30+0.03 47.413. 28.414. 72.431. -0.35+40.42 80,7 IRS 13E-WP? 
HCRSI 39.7672(1) 31.378(3) 0.19+0.01 0.24+0.02 unresolved —0.78+ 0.40 7.0, 7 IRS 2-W 
HCRS2 39.7620(4) 29.783(4) 0.3140.04 0.4740.03 24.41. 1641. 60.43. -1.38+0.49 8.0, 7 IRS 13E-W!? 
HCR53 39.7513(1) 30.675(3) 0.36+0.02 0.3040.03 80.411. 14.441. 20.43. 0.6140.34 7.0, 9 IRS 13E-SW”? 
HCR54 39.7461(3) 26.567(10) 0.1140.02 0.15+0.02 28.410. 12.421. 60.410. -1.03+40.78 7.0, 8 IRS 34SW” 
HCR55 39.7439(4) 31.215(6) 0.5740.08 0.86+0.14 81.422. 53.446. 108.498. -1.3740.69 7.0, 7 
HCR56 39.7320(1) 30.580(4) 0.19+0.02 0.3640.02 58.47. 38.46. 55.410. -2.1240.43 7.0, 8 
HCR57 39.7299(1) 26.595(3) 0.33+0.03 0.35+0.03 38.417. 18.426. 60.442. -0.20+0.48 7.0, 8 IRS 34SW~”, ID6”> 
HCR58 39.6399(1) 26.651(4) 0.19+0.02 0.20--0.02 76.418. 44.446. 124.453. -0.1740.54 7.0, 8 IRS 6E? 
HCR59 39.6255(3) 26.623(4) 1.2140.08 1.404017 94.422. 68.47. 132.435. -0.48+0.47 80,8 IRS 6E”, ID37?5 
HCR60 39.5463(1) 35.005(1) 0.60+0.04 0.50+0.03 31.412. 16.450. 90.423. 0.6140.30 7.0, 8 D75 
HCR61 39.4593(1) 31.736(4) 0.34+0.03 0.31: 0.03 unresolved 0.3140.43 7.0, 9 ID85 
HCR62 39.3096(1) 30.683(1) 0.22+0.01 0.18+0.01 30.45. 9.42. 10.42. 0.672033 7.0, 8 
HCR63 392771(4) 27.332(12) 0.22+0.04 0.27+0.02 51.433. 36.434. 42.448. -0.68+0.69 7.0, 8 
HCR64 40.0201(9) 31.178(20) <0.15 0.60+0.07 140.420. 70.410. 134.+10. NA 8.0, 5 micro-qso,K4074 
1 The radio counterpart of SGR J1745-2900 |Eatough et al. [2013], the GC magnetar, and see Sec.4.3 of this 2 peated east of the IRS 7 bow-shock feature [Yuset-Zadeh & Moris [ooi] . 3A candidate cannonball of the PWN G359.95-0.04 
Wang et al. [2006], and see Sec.4.1 of this paper. ^ The head of the bullet [Zhao et al. | onents. 2 The radio core associated with IRS 13E. Located ~0.5" SE of IRS 13E core (e.g. | Tsuboi et al. 2017|. "Lo 


Te: 
the mini-cavity. "The e source that has been resolved into three components | Yusef-Zadeh et al. [1 1991], S (ŒHCR17), NE and NW CHCR20) the three components correspond to RS6, RSS and RST detected at 34.5 GHz 
6]. Loca ated in the northern mini-cavity. 101 oca ated at the northwestern rim of the mini-cavity. ^ Located 0.5" NE of the IRS 13E radio core. 2 Located at the southwestern rim of the mini-cavity. 3 Located 1" NE of the IRS 13E radio 
core. Located at the western rim of the mini-cavity. !5 Located 0.7" N of the IRS 13E radio core. !®Located N of the IRS 2. !7 Located 0.2" S of the IRS 13E radio core. 19A triplet located 1.2" N of the IRS 13E radio core, consisting of three 
components S(HCR40), NW(HCR43) and NE with a size larger than the HCR upper limit. !? Located 0.3" W of the IRS 13E radio core, the HCR41 group consisting of six members of HCR41, 44, 46, 48, 50 and 52. 20 Located 1" SW of the IRS 13E 


and see Sec.4.2 of this paper. 25 The 7mm-sources are identified with IR stars emitting at 3.8 pm, among which HCROS (ID4), HCR06 (ID2), HCR10 (ID1), HCR57 (ID6), HCR60 (ID7) and HCR61 (ID8) are associated with strong stellar 


Yusef-Zadeh et al. i 


winds 


4.3. The GC magnetar: SGR J1745-2900 


HCR07 is one of the most recognizable HCR members, as 
it is associated with the GC magnetar and has a high-energy 
counterpart, SGR J1745-2900. The soft gamma-ray repeater 
was discovered by Swift during a large X-ray outburst on 2013 
April 24 (MJD 56406), powered by a magnetar close to Sgr 
A*(Kennea et al. [2013). The magnetar hypothesis was further 
supported by NuSTAR detections of periodic pulsed signal at 
3.76 s Tees Jpory, With the observations by Chandra and 
Swift, (2013) pinpointed the location of the magnetar 
at a projected distance of 2.443-0.3 arcsec from Sgr A*; and the 
authors also determined the source spin period and its derivative 
with high precision (P 23.7635537(2) s and P= 6.61(4) x 107? 
s s !). The magnetar, SGR J1745-2900, was monitored by the 
Chandra X-ray observatory for six years following the X-ray 
outburst in 2013 April, showing the long-term properties of the 
outburst (Rea et al. [2020). 

Radio pulses from SGR 71745-2900 were first detected by 
the Effelsberg radio telescope and confirmed with other tele- 
scopes, including the Nancay telescope, JVLA, Jodrell Bank 
Observatory TATUR) Sho RE Topmsion POT at varichs 
Compact Array (ATCA) (Shannon & Johnston |2013) at various 
frequencies between 1.5 and 19 GHz. High-frequency pulses 
were detected at 87, 101, 138, 154, 209 and 225 GHz with 
the IRAM-30m telescope during the period between 2014-7- 
21 (MJD 56859) and 2014-7-24 (MJD 56862) (Torne et al. | 
(2015). In their follow-up campaign, 
tected high-frequency pulses with the IRAM-30m up to 291 
GHz in the interval between 2015-3-4 (MJD 57085) and 2015- 
3-9 (MJD 57090). The high-frequency pulses were also de- 
tected at 45 GHz with the Green Bank Telescope (GBT) on 


2014-4-10 (MJD 56757) (Gelfand et al. |2017). 


4.3.1. A collection of flux-densitiy measurements from JVLA 
and ALMA 


In addition to the JVLA and ALMA flux densities deter- 
mined from the images reported in this work, we also collected 
the data from the prior published literature. Table 3 assembles 
all the available data from JVLA and ALMA observations at ra- 
dio wavelengths. The table is configured into two main column 
sections, each of which consists of six sub-columns: (1) observ- 
ing date; (2) the corresponding modified Julian day (MJD); (3) 
band center frequency in units of GHz; (4) flux densities in units 
of mJy and (5) the corresponding 1c uncertainties. For the non- 
detections, a 3c value is given for the upper limits. Finally, al- 
phabetical codes are designated for the relevant references that 
are listed at the bottom notes of the table. The measurements 
reported for the first time in this paper are highlighted with a 
bold font. 


4.3.2. Radio spectrum of the GC magnetar 


In addition to the JVLA and ALMA data, we also collected 


the flux-density data of the magnetar from ATCA (Shannon &| 
[Johnston |2013), yielding a total of 91 flux-density measure- 


ments of this object at radio wavelengths over the 6.5 yrs since 
the 2013 outburst. To investigate the spectrum of the magne- 
tar, we binned the data into nine bands with frequency ranges 
corresponding to those used for the ALMA and JVLA observ- 
ing bands. We then computed the mean flux density (S,) and 
variance (o2) in each of the bands, weighted by wt(i) = gor? 
where o(i) is the uncertainty of each individual flux density, i. 
For non-detections, a zero weight was adopted, assuming that 
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Data: JVLA, ALMA & ATCA 
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FIG. 6.— Averaged spectrum of SGR J1745-2900, made with band-averaged 
flux densities observed with the JVLA and ALMA (values summarized in Ta- 


ble 3) as well as with ATCA (Shannon & Johnston |2013). 


the actual uncertainties of the non-detections due to the errors in 
the calibration for the system and atmostheric issues are much 
greater than the cited rms errors. For a total of n measurements 
in each bin, the number of non-zero weighted data points, m, is 
less than or equal to n, (m <n), within each bin. The error of 
the mean flux density can be determined with 


2; 
Hs mA 

We note that the Ka band shows a lower detection rate (~55%) 
than the Q-band (768406), indicating that a spectral minimum is 
located near the Ka-band frequency of 33 GHz, given the fact 
that the rms noise at the Ka-band center frequency of 33 GHz is 
about a factor of 2 smaller than that of Q-band (see Table 1). To 
reduce the bias due to zero-weighting on the data with negative 
detections, we combined the Ka-band and Q-band bins’ data 
and re-computed the mean flux density and its uncertainty. Ta- 
ble 4 summarizes the results. Fig. 6 shows the averaged spectra 
observed with the JVLA, ALMA and ATCA over 6.5 yrs since 
the outburst. The vertical bars on each point reflect the large 
intrinsic variation of the flux density from SGR J1745-2900 
over the observed time scale. However, a spectral minimum 
appears to be present at a frequency around 30 GHz, hereafter 
referred to as the transition frequency , which appears to sep- 
arate the spectrum into two components arising from two emis- 
sion regimes at centimeter and millimeter-submillimeter wave- 
lengths. The flux density of the centimeter-wave component 
is typically ~1 mJy while the millimeter-wave component is 
about three times more intense. 

We then carried out a least-squares fitting of the spectrum 
to both the cm- and mm-components with two power-law func- 
tions. The centimeter data are well fit with a power-law spectral 
index of aq =—1.5+0.6 (S, x vy), a steep spectrum similar 
to that of radio pulsars. A steep spectrum of o — —1 was pre- 
viously reported between 4.5 and 8.5 GHz for the observations 


during the 2013-5-30 flare of the cm-component (Shannon &| 
[Johnston [2013). The spectral data points rising through mil- 


limeter wavelengths can be fit with a spectral index of oq = 
1.1 x: 0.2, indicating the presence of an emission bump or a 
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TABLE 3 
JVLA & ALMA FLUX DENSITIES OF SGR J1745-2900 


Obs. date MID  vo(GHz) S, (mJy) o(mJy) Ref. 


Obs. date MJD vo (GHz) S, (mJy)  e(mJy) Ref. 


a) 2 (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 
2011-08-04 55776 42 « 0.24 a 2012-10-14 56214 212 « 0.60 a 
2012-10-14 56214 32 < 0.81 a 2012-10-14 56214 41 < 0.78 a 
2012-12-22 56283 21.2 < 0.90 a 2012-12-22 56283 32 < 0.94 -— a 
2012-12-22 56283 41 < 0.69 — a 2013-05-10 56422 9 0.56 0.011 g 
2013-06-01 56444 9 0.76 0.015 g 2013-06-30 56473 15 0.58 0.012 g 
2013-07-13 56486 9 1.47 0.030 g 2013-10-26 — 56591 21.2 < 1.95! a 
2013-10-26 56591 32 < 1.45 isa a 2013-10-26 56591 41 < 0.82! a 
2013-10-26 56591 41 0.7 0.4 e 2013-11-29 56626 21.2 < 2.04" a 
2013-11-29 56626 32 < 0.82 -- a 2013-11-29 56626 41 < 0.70! a 
2013-11-29 56626 41 0.89 0.08 e 2013-12-29 56656 21.2 < 1.707 a 
2013-12-29 56656 32 < 0.86 = a 2013-12-29 56656 41 < 1.52! pee a 
2013-12-29 56656 41 1.20 0.7 e 2014-01-01 56658 15 2.09 0.042 g 
2014-01-01 56658 9 1.18 0.024 g 2014-02-15 — 56703 21.2 0.84 0.33 a 
2014-02-15 — 56703 32 1.83 0.10 a 2014-02-15 — 56703 41 .85 0.07 a 
2014-02-15 56703 41 2.1 0.4 e 2014-02-21 56709 44.6 .62 0.04 a 
2014-02-22 56710 15 1.07 0.022 g 2014-02-22 56710 9 0.94 0.019 g 
2014-03-09 56731 34.5 1.30 0.01 a 2014-03-22 56738 21.2 2.79 0.19 a 
2014-03-22 56738 32 2.64 0.05 a 2014-03-22 56738 41 .24 0.02 a 
2014-03-22 56738 4l 2.1 0.3 e 2014-04-03 56750 23 0.92 0.019 g 
2014-04-03 56750 43 0.54 0.011 g 2014-04-25 56772 9 .00 0.020 g 
2014-04-25 56772 15 1.22 0.025 g 2014-04-26 56743 21.2 0.90 0.14 a 
2014-04-26 56743 32 0.62 0.04 a 2014-04-26 56743 41 .20 0.07 a 
2014-04-26 56743 41 0.91 0.30 e 2014-04-17 56764 9.0 3.50 0.08 b 
2014-05-10 — 56787 41 1.15 0.05 e 2014-05-17 56794 5.5 4.50 0.24 b 
2014-05-26 56803 5.5 3.90 0.09 b 2014-05-31 56808 21.2 4.21 0.17 a 
2014-05-31 56808 32 2.90 0.13 a 2014-05-31 — 56808 41 2.94 0.12 a 
2014-05-31 56808 41 3.5 0.4 e 2014-08-23 56892 15 0.63 0.013 g 
2014-08-30 56899 23 0.26 0.006 g 2014-08-30 56899 43 0.15 0.003 g 
2015-02-20 57073 9.0 3.00 0.30 b 2015-09-11 57276 33.0 1.80 0.05 b 
2015-09-16 57281 44.6 2.45 0.05 b 2016-04-23 57502 343 2.80 0.23 b, d 
2016-07-12 57581 44.2 5.79 0.05 f 2016-07-15 57584 226 4.70 1.29 f 
2016-08-30 57630 343 3.10 0.11 b.d 2016-08-31 57631 343 3.29 0.12 b.d 
2016-09-08 | 57639 343 3.90 0.13 b.d 2017-10-06 — 58032 225 1.41 0.08 b.c 
2017-10-07 — 58033 225 1.53 0.08 b.c 2017-10-09 — 58035 225 0.79 0.07 b.c 
2017-10-10 — 58036 225 1.22 0.06 b.c 2017-10-11 — 58037 225 1.25 0.07 b.c 
2017-10-12 58038 225 1.44 0.07 b.c 2017-10-14 58040 225 3.16 0.08 b.c 
2017-10-17 58043 225 5.34 0.09 b.c 2017-10-18 — 58044 225 5.42 0.11 b.c 
2017-10-20 58046 225 5.60 0.12 b,c 2019-09-08 58734 5.5 < 0.507 sigs b 
2019-09-21 58747 9.0 0.40 0.10 b 2019-09-20 58748 318 1.32 0.15 b 
2020-11-27 59180 9 0.36 0.06 b 


Detected in 


TABLE 4 


MEAN FLUX DENSITIES OF SGR J1745-2900 IN THE 
JVLA/ALMA BANDS 


Band code n m VXAv(GHz2 S 7g (mJy) 


a) 2 (3) (4) (5) 
A7 5 5 338.0+5.0 
A6 11] 11 225.1+ 0.1 
Q-Ka 34 23 39.4+ 0.9 
K 15 10 20.7+ 0.5 
Ku 7 7 15.3+ 0.2 
X 12 12 8.9+ 0.1 
C 7 6 5.5+ 0.4 


The band-averaged data are derived from the individual flux-density mea- 
surements with the JVLA and ALMA (Table 3) as well as ATCA 
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plateau at higher frequencies. Such a high-frequency plateau 
appears not to be unique to the GC magnetar. The radio-active 
magnetar 1E 1547.0—5408 has also been observed to have a 


.0 
T A 3-a upper limit. Reference: a. |Yusef-Zadeh et al. (2015). b. This paper. c. The ALMA observations at the ten epochs were carried out by the PI Tsuboi 
in (Tsuboi et al. |2017). We re-processed the ALMA archival data (2017.1.00500.S) and determined the flux densities of SGR 11745-2900 at 225.75 GHz. d. 
e ALMA image made from observations on 2016-4-23, 2016-8-30/31, and 2016-9-8 (Tsuboi et al. |2017). We re-processed the ALMA archival 
data (2015.1.01080.S) and determined the flux densities of SGR J1745-2900 from each of the 340-GHz images a 


from two VLA Q-band measurements with a baseline filter (7 500 kA) and all data 
measurements is adopted. f. |Yusef-Zadeh et al. (2017). g. 2015/16). 


€ four epochs. e. Mean values determined 


Gelfand et al. |2017); for the case of std = 0, the smaller error in the two 


spectrum rising towards millimeter wavelengths (Chu et al. | 
2021). Overall, the spectrum of the GC magnetar at the fre- 


quencies in the range between 5 and 310 GHz is described with 
a combination of the two power-law functions (see Fig. 6) with 
a transition frequency v; ~ 30 GHz corresponding to a mini- 
mum flux density of ~0.3 mJy. 


4.3.3. Radio variability of SGR J1745-2900 


We inspected the radio variability of the GC magnetar SGR 
J1745-2900, including a total of 161 flux-density measurements 
since the onset of the 2013 outburst. The measurements by 


(2015| [2017) are also included in addition to 


the JVLA-ALMA (Table 3) and ATCA (Shannon & Johnston 
[2013) data. The magnetar appears to be highly variable on all 


observed time scales and wavelengths (e.g. 
and this paper). We uti- 
lized a bin-averaging algorithm similar to the analysis of the 
magnetar spectrum (Sec.4.3.2). We binned the MJD or the time 
axis with a constant interval of 100 days. The two spectral com- 
ponents —— low-frequency and high-frequency — correspond- 
ing to the cm-component (v < 30 GHz) and the mm-component 
(v > 30 GHz), were examined separately. We then computed 
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FIG. 7.— (a) The X-ray luminosity variation after the onset of outburst 


of SGR J1745-2900, observations by the Chandra X-ray observatory 
al. |2020). (b) 100-MJD bin averaged radio light curve of the magnetar 

900 for both the mm- (red dots) and cm- (black dots) components. In addition 
to the flux-density measurements from the JVLA and ALMA summarized in 


Table 3, other telescope data from literatures ee NC SUE 
[Tome eral 2015] 2017) are included. (c) Detail o MA observations 
in the interval between 2017-10-6 and 2017-10-20 (1626 days after the out- 
burst of SGR 1745-2900 on 2013-4-23) show a day-to-day variability of the 
magnetar SGR J1745-2900 at 225 GHz displayed as normalized flux densities 
SMID-56406/Smin (red dots), where the minimum flux density Smin =0.79+0.08 
mJy for the magnetar on MJD 58035. The variability of SGR J1745-2900 dur- 
ing the two-week interval is compared with that of Sgr A* (dark-blue dots); 
Sgr A* (Smin=3.015+0.013 Jy and Save= 3.26+0.09 Jy) that is three orders in 
magnitude more luminous than the magnetar. 


the weighted mean flux densities and the corresponding dis- 
persions in each of the 56 MJD-bins for both mm- and cm- 
components. The results for the non-empty bins are tabulated 
in Table 5. Columns 1 and 2 show the bin ID (binID) and the 
corresponding central MJD. Columns 3 and 4 are the numbers 
of all the measurements (n) and the non-zero weight data (m) 
in the corresponding MJD-bins; as above, the measurements 
with only upper limits are zero-weighted. Column 5 is the 
mean MJD of the observing dates. Column 6 gives the mean 
flux-density (SbintD)s weighted by the inverse variance on each 
of the measurements, along with the uncertainty of the mean 


V Sub 
dispersion due mainly to the variation in flux density, is listed 
in column 7. 

Fig. 7b plots the MJD-bin averaged radio light curves for the 
two spectral components, red for mm and black for cm. Unlike 
the X-ray light curve (Fig. 7a) observed by the Chandra X-ray 


=\ on /m. And the standard deviation (pinn), or the 
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observatory (Rea et al. |2020) that shows a smooth decrease 


in X-ray luminosity, both the cm- and mm-components varied 
significantly in the first 800 and 1700 days. Our two measure- 
ments of flux density based on the earlier JVLA observations 
at 5.5 GHz showed a significant variation from 4.50+0.24 mJy 
on 2014-05-17 to 3.90+0.09 mJy on 2014-05-26, on the time 
scale of a week, which is consistent with the large variability at 
5 GHz observed during May 2013 (Shannon & Johnston 2013). 

The two most recent JVLA measurements at 9 GHz on 2020- 
11-27 (MJD 59180) and 2019-09-21 (MJD 58777) indicate that 
the flux density of the cm-component decreased to 0.4 mJy 
from 3.5 mJy in the JVLA observation on 2014-04-17 (MJD 
56764). Also, we failed to detect the source at 5.5 GHz based 
on the JVLA observations on 2019-09-08 (58734), imposing 
a 3-c upper limit of 0.5 mJy while the flux density was 4.5 
mJy based on our JVLA observation at 5.5 GHz on 2014-5-17 
(MJD 56794). The mean flux density of the cm-component of 
the magnetar decreased by a factor of 3 over a period of 6.5 yr 
(Fig. 7b and Table 5). 

The Chandra X-ray observations of SGR J1745-2900 follow- 
ing the outburst show a smooth decrease from 4.9x 10? erg ^! 
at the onset of the X-ray outburst to 0.047 x 10?? erg ^! at the 
most recent epoch (Fig. 7a). The X-ray luminosity therefore 
dropped by two orders in magnitude over 6 years. The general 
trend of decreasing radio flux density is consistent with the X- 
ray light curve, although the cm-component has a much slower 
decline, is much more variable, and shows a large range of vari- 
ability on time scales from days to years. 

The mm-component stayed at a relatively low level, with bin- 
averaged flux densities of 0.7+0.40 and 1.53+0.23 mJy dur- 
ing the second and third MJD-bins after the outburst (Table 5 
and Fig. 7b). During the first 300 days, the mm-component 
was difficult to detect (only 5 detections out of 11 observa- 
tions). We note that the JVLA Q-band data observed during 
2013 fall and 2014 spring were reduced by two independent 
groups (Yusef Zadeh e. DUIS] Gelfand et al. 2017. Th 
mm- E reached a maximum of 5.63+0.19 mJy in the 
7th MJD-bin (601-700 days after the outburst), based on 35 ob- 
servations with a detection rate of 100%. 

The 17th MJD bin of the mm-component (between 1600- 
1700 days) contains 10 individual ALMA observations at 225 
GHz, all carried out within a two-week period during October 
2017 with a high-resolution (0.02") configuration 
(2019). The mean flux density for the mm-component in this bin 
was 2.25+0.56 mJy. With this configuration, the extended HII 
emission, as well as dust emission from the local medium, are 
well resolved out and only hyper-compact sources can be de- 
tected. The typical rms noise of the ALMA images is about 
20-30 uJy beam !, and the magnetar has S/N ratios of 50-100. 
At a distance of ^3" from the magnetar, Sgr A*, with a flux 
density near 3 Jy, is an excellent reference source to examine 
the day-to-day variability of the magnetar. The ALMA data 
show the magnetar to be highly variable at 225 GHz on a time 
scale of days. The source started at 1.41-2-0.08 mJy on 2017- 
10-6, dropped to a minimum of 0.79+0.08 mJy three days later 
on 2017-10-9, and then increased by a factor 7 within 8 days 
reaching 5.34+0.09 mJy on 2017-10-17, staying at that level 
for the next several days. To compare the variability of SGR 
J1745-2900 with Sgr A*, we normalized the source flux densi- 
ties by their own minimum flux densities during the observing 
period. Fig. 7c shows the relative variability for both the mag- 
netar and Sgr A* during the two weeks in bin 17. We define 
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a relative variability parameter, ZY, to quantitatively describe 
the magnitude of variability relative to a minimum flux density 
Smin, given a maximum flux density Smax of a target source ob- 
served in a period: 


[Smax x Smin] 

Simin ` 
The magnetar shows the relative variability ZV ~ 6 while Sgr 
A* is observed to be only moderately variable in the same ALMA 
program, with ZY ~ 0.2.. 

Finally, the latest available detections of the magnetar in 2019- 
2020, are 1.32+0.15 mJy at 318 GHz with ALMA on 2019-09- 
20 (MJD 58748) and 0.36+0.06 mJy with the JVLA on 2020- 
11-27 (MJD 59180). 


BV = 


TABLE 5 


MEAN FLUX DENSITIES OF SGR J1745-2900 IN 100-DAY MJD 
BINS 


binID MJD nm MJD 
(day) (day) 


a) 2 6) 4 (5) (6) (7) 


Soup £ 0: TbinID 
Spin 


(mJy) (mJy) 


Centimeter-wave component 


1 56456 16 16 56435.1 0.50+0.06 0.24 
3 56656 5 3 56673.2 1.3640.26 0.45 
4 56756 10 10 56755.7 1.08+0.14 0.44 
5 56856 12 12 568624 0.46+0.09 0.31 
7 57056 10 10 57085.9 0.79+0.13 0.41 
24 58756 2 1 58777.0 0.40+0.10 0.10 
28 59156 1 1 59180.0 0.36+0.06 0.06 
Millimeter-wave component 
2 56556 3 1 56591.5 0.70+0.40 0.40 
3 56656 9 5 56678.6 1.5340.23 0.51 
4 56756 10 10 567424 1.00+0.14 0.44 
5 56856 26 21 568532 2.203021 0.96 
7 570566 37 37 570873 5.6340.19 1.16 
9 57256 2 2 57278.5 2.1340.32 0.45 
11 57456 1-1 57502.0 2.80+0.23 0.23 
12 57556 2 2 57582.5 5.7840.04 0.06 
13 57656 3 3 57633.3 3.40+0.23 0.40 
17 58056 10 10 58038.4 2.14+0.53 1.68 
24 58756 1 1 58770.0 1.32+0.15 0.15 


The 100-day-averaged data are derived from a total of 161 individual flux-density 
measurements with the JVLA, ALMA (Table 3), ATCA (Shannon & Johnston 12013 


and other single-dish telescopes 20151|2017] since the onset of the 
outburst at the epoch 2013-4-23 (MJD 56406). 


5. ASTROPHYSICAL IMPLICATIONS 


A population of hyper-compact radio sources (HCRs) are de- 
tected at 33 and 44.6 GHz with the JVLA in the vicinity of Sgr 
A* within a radius of 13 arcsec. The new survey was motivated 
by the previous JVLA detections of 110 Galactic center com- 
pact radio sources (GCCRs) at 5.5 GHz in the radio bright zone 
within a radius of 7.5 arcmin from Sgr A* but outside Sgr A 
West. 


5.1. Spectral types of HCRs & their distribution in 
flux-density 
Fig. 8 shows the flux-density distribution of the HCRs as 


compared with that of the GCCRs. The distribution of the GC- 
CRs is similar to the high-luminosity tail of the pulsars’ distri- 


bution in the Galactic disk (Kramer et al. [1998 


Zhao, Morris & Goss |2020). The HCRs have a rel- 
atively narrow distribution from —1 to 0.6 in Log(S [mJy]), or 
ranging from 0.1 to 4 mJy, peaked at —0.4 in Log(S [mJy]) or 
0.4 mJy. We note that the flux density of the peak in the HCR 
distribution is close to the minimum value 0.32 mJy at the tran- 
sition frequency observed in the band-averaged spectrum for 
the GC magnetar J1745-2900 (see Sec.4.3.2). 

We divided the HCRs into three sub-types according to their 
spectral indices a44.6/33 (S x v^**5/5) derived between 44.6 and 
33 GHz: flat (0.2 > 044 6/33 > —0.5), steep (044.5/33 < —0.5) and 
inverted (044,533 > 0.2). 

Of all the HCRs, 58% (38/65) are steep-spectrum sources, 
26% (17/65) are inverted-spectrum sources, and 15% (10/65) 
have a flat spectrum. We note that HCR49 is a double; so a total 
of 65 HCR components are included in the spectral-index anal- 
ysis. The inset of Fig. 8 shows the flux-density distributions for 
each of the three spectral types with a finer bin, ALog(S [mJy]) 
- 0. 


5.1.1. Flat-spectrum HCRs 


The flux density distribution of the flat-spectrum HCRs ap- 
pears to be uniform within the range from 0.16 to 1.6 mJy. 
Some fraction of the flat-spectrum HCRs might be unresolved 
peaks in the HII region components since a flat spectrum at Ka- 
and Q-bands is characteristic of free-free emission in optically 
thin HII regions. 


5.1.2. Steep-spectrum HCRs 


This sub-sample consists of 38 members, the largest sub- 
sample among the three, in which the flux densities are sta- 
tistically well distributed. The distribution can be fitted with a 
Gaussian, with a mean value of u = —0.35 and a standard de- 
viation of o = 0.22 in Log(S [mJy]). The mean value of the 
steep-spectrum HCRs corresponds to 0.45 mJy in flux density. 

The spectral index 0446/33 of the steep- spectrum HCRs is in 
the range between —6.5 and —0.6, giving a mean Ha =—1.8+ 
0.2. The steep spectrum of this sub-sample differs distinctly 
from the HII components, suggesting the presence of a popula- 
tion of hyper-compact nonthermal radio sources in the central 
parsec. The nonthermal HCRs are likely associated with the 
massive stellar remnants that are expected to be distributed in 


the close vicinity of Sgr A* (Morris |1993}|Hailey et al. |2018} 
Generozov et al. 2018). 


5.1.3. Inverted-spectrum HCRs 


For this sub-sample, the mean flux density and standard de- 
viation are u = 0.5 mJy and ø = 0.25 mJy, respectively. The 
spectral index 0445/33 of the inverted-spectrum HCRs is in the 
range between 0.21 and 1.65, with a mean value of pa 2 0.61 + 
0.08. Among the three spectral sub-samples, the distribution of 
the inverted-spectrum HCRs appears to most closely follow the 
GCCR distribution, which matches the high-luminosity tail of 
the distribution that normal pulsars would have at the Galactic 
center (Zhao, Morris & Goss |2020). Normal pulsars usually 
have a steep spectrum at centimeter wavelengths and are diffi- 
cult to detect at high frequencies. However, the GC magnetar, 
emitting at the Ka- and Q-bands, falls into this sub-sample. As 
shown in Sec. 4.3.2, an inverted spectral component of SGR 
J1745-2900 is present at millimeter wavelengths in the aver- 
aged spectrum of a large sample of observations. The inverted 
spectrum cannot be simply attributed to time variability of the 
flux density. By analogy with SGR J1745-2900, the inverted 


spectral component could be indicative of an association of the 
HCRs with magnetars. 

On the other hand, hyper-compact HII regions associated 
with a late-type massive star may also have a similar character- 
istics in spectral index, flux density level and compactness. For 
example, a group of compact components have been detected at 
mm-submm wavelengths in the IRS 21 complex (see Fig. 2b). 
The hyper-compact sources associated with IR emission more 
likely belong to hyper-compact HII or HC HII regions associ- 
ated with young stellar objects rather than old massive stellar 
remnants. The hyper-compact radio components of IRS 21 are 
therefore not included in the HCR catalog (Table 2) in this pa- 
per. IRS 21 will be discussed in a separate paper. 

From cross-correlation between the HCRs of this paper and 
the nine 7mm-IR sources that are found to be associated with 
strong stellar winds (Yusef-Zadeh et al. |2014), we find that 
six of the nine candidate IR stars have HCR counterparts. The 
spectral indices of the six HCRs - five with inverted spectra and 
one with a flat spectrum — are consistent with the radio emission 


produced by the ionized winds of hot, massive stars (Panagia &| 


1975). Thus, a large fraction of the 17 inverted-spectrum 
HCRs might consist of thermal free-free emission sources. Ap- 


parently, a significant portion of the inverted-spectrum HCRs is 
associated with late-type massive stars. 

We speculate that more magnetars besides SGR J1745-2900 
reside in our inverted-spectrum sample of HCRs. However, 
we must be able to distinguish such objects from compact ra- 
dio sources associated with young massive stars. Further high- 
resolution ALMA observations of the variability and polariza- 
tion characteristics of this subsample will be crucial for identi- 
fying the nonthermal nature of candidate magnetars and stellar- 
mass black holes. 


5.2. A dense group of HCRs and the radial distribution of 
their surface-density 


The surface density of 110 GCCRs, located outside Sgr A 
West (> 1 arcmin from Sgr A*) but within the RBZ (< 7.5 ar- 
cmin from Sgr A*), is ~0.6 counts arcmin ?. This surface den- 
sity is an order of magnitude greater than that of background 
extragalactic compact radio sources (e.g. 
Gim et al. |2019). The dense group of 64 HCRs, located within 
a radius of 13" from SgrA*, has a relatively higher surface den- 
sity, with an average value of ~500 counts arcmin ?. We con- 
clude that contamination of our sample by extragalactic back- 
ground sources is very likely to be negligible. 

In Appendix A, we developed a procedure to construct the 
surface-density (Xers) distribution of the radio compact sources 
(RCS) as function of the projected radial distance (R) from Sgr 
A*. including 64 HCRs at R « 13" (this paper) and the 110 GC- 
CRs outside Sgr A West (Zhao, Morris & Goss |2020) as well 
as the Galactic center transient (GCT, see [Zhao et al. |1992). 
Fig. 9 shows the radial distribution of the surface-density of 
the RCSs detected within the RBZ. Excluding the four data 
points related to a Xers—bump in the radial distance range be- 
tween 1.5" and 7", the surface-density distribution can be fitted 
with two power laws (see Fig. 8) as described by the “Nuker” 
law (Lauer et al. [1995 
2018). In the inner region at R « 30" but excluding 
the range between 1.5" and 7", the surface-density distribution 
of the HCRs follows a power law: 


Eycr(R <30”) ~ (5.02 1.8) x 10°R 2... (arcmin) (1) 


arcsec 


13 


Histograms of HCRs—33GHz vs GCCRs—5.5GHz 
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FIG. 8.— Flux-density distribution of the HCRs at 33 GHz within the central 
13" x 13" region versus the GCCRs detected in 
the GC radio bright zone (15' x 15"). The grey histogram represents the 110 
GCCRs detected at 5.5 GHz outside Sgr A West. The 64 HCRs (65 compo- 
nents, HCR49 is double), shown by the hatched histograms, are divided into 
three types: flat —0.5 < a < 0.2 (red), steep a < —0.5 (green) and inverted 
a > 0.2 (blue), according to their radio spectral index determined between 
44.6 and 33 GHz, except for three steep-spectrum sources - HCR22, HCR32 
and HCR64 — which were not detected in the Q-band data used in this paper. 
The spectral indices for HCR22 and HCR32 are determined using X-band data 
discussed in this paper. The spectral index for HCR64, the micro-quasar, is de- 
termined using the K-band data p TED uus assuming that the source 
was in a quiescent state during the 2015 observations. The inset panels on 
the right are the more finely binned individual distributions of the three spec- 
tral sub-samples. The distribution of the steep-spectrum HCRs is fitted with a 
Gaussian (red curve). 
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FIG. 9.— Radial distribution of the surface density of compact radio compact 
sources detected at the Galactic center. The red points are the surface density 
data X% corrected for a higher equivalent sensitivity cutoff at 5.5 GHz by 
multiplying the original 33-GHz surface density data, Uers (light-grey points), 
by a correction factor, 7 = 2.4, derived from the 33-GHz observations of the 
HCRs (see Appendix A of this paper). The black points mark the data of 
GCCRs observed at 5.5 GHz. In the annular zones around R = 3”, a bump in 
Xs is present. The orange line is the best fit using a power-law with an index 
T = 1.6 30.2 for the inner region outside of the bump at 1.5 - 7 arcsec; and 
the blue line indicates the least squares fitting to the steeper power law with 
B - 2.0 0.2 for the outer region. The break radius, R}, is ~ 30", or ~ 1.2 
pc. The top-right inset shows the surface-density distribution of the HCRs 
excluding the flat-spectrum sources. 
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where the index [= 1.6+0.2. At large radii R > 30", the 
surface-density distribution of the GCCRs shows a steeper power 
law: 


(arcmin) (2) 


where 8 =2.0+0.1. At R ~ 30", the two power laws intersect. 


Xacca(R > 30”) ~ (2.0 4 0.9) x 104 R78 


arcsec 


We note that the power-law distribution derived from HCRs 
is much steeper than that I ~ 0.6 of the stellar cusp determined 
with the K, stars in magnitude range between 12.5 and 18.5 in 
the radial distance range between 1" and 50" 


al. [2018). 


5.3. Candidate massive stellar remnants 
5.3.1. Mass segregation 


It has been known for several decades that the central parsec 
contains a large number of young massive stars (Krabbe et al. | 
1995} 
, with ~100 
O and Wolf-Rayet (WR) stars confined within a radius of ~0.4 
pe Ghez et al. [2005). These stars are rela- 
tively young (< 10 Myr) and are orbiting Sgr A*, the supermas- 
sive black hole (SMBH) at the Galactic center. The early-type 
stars provide substantial UV photons to maintain the ionization 
of the gas within the central parsec (Zhao et al. [2010). 

The stars in the central cusp apparently have a top-heavy 
present-day mass function (over-abundance of high-mass stars, 
see [Lu et al. 2013); and the fact that the 
orbits of many of them collectively define a coherent disk sug- 
gests that the young nuclear cluster may have originated via in 


situ star formation (Levin & Beloborodov |2003 2007 
Nayakshin et al. [2007,|Collin & Zahn |2008). The process that 


formed the young nuclear cluster near the central black hole is 
likely to have been a recurring one. This process follows a limit 
cycle of activity wherein the star formation is a violent event co- 
inciding with a heavy accretion episode onto the SMBH. Thus 
the inner disk is quickly disrupted and dissipated in the process. 
Following the disruption of the inner disk, continued migration 
of gas from the central molecular zone toward the center re- 
builds the disk, eventually leading to the next star formation 
event (Morris, Ghez, & Becklin |1999). With repeated instances 
of this cycle, the remnants of massive stars are produced and 
left in place to collect at the bottom of the Galaxy's gravita- 
tional potential well. Furthermore, the more massive remnants, 
particularly stellar-mass black holes, will undergo mass segre- 
gation in the dense central star cluster, as a result of dynamical 
friction, and become even more concentrated toward the black 


. The mass ere ene may partially explain the 
RENT flux of gamma rays and X-rays from the Galactic center 
caused by the possible presence of a large population of mil- 
lisecond pulsars (MSPs) and cataclysmic variables (CVs), as 
suggested by (Arca-Sedda et al. [2018). The excess of GeV 
gamma rays toward the Galactic center may alternatively be ex- 
plained by a high supernovae rate, leading to the production of 


neutron stars and ultimately to a MSP population (O'Leary et| 
al. [20152016] |Calore et al. |2016). 


In the case of stellar-mass BHs that are significantly more 
massive than the mean stellar mass expected for an evolved 


population, such heavier objects would migrate toward the cen- 
ter and be distributed in a compact cluster around the SMBH. 
The nature of the mass segregation depends on the relaxational 
coupling parameter A = 4ngyMĝy/ [ng M2 G * Msu/ M)] for 
black hole (BH) mass Mgy and spatial number density ngu 
along with stellar (ST) mass M, and spatial number density 
ng. The value of A is a measure of the importance of BH-BH 
scattering relative to BH-ST scattering for the dynamical fric- 
tion process (Alexander & Hopman [2009). Weak segregation 
occurs when A >> 1; and when A < 1, the strong segrega- 
tion applies. The latter case leads to a steeper slope in the 3D 
radial density distributions of BHs and a larger central concen- 
tration of BHs relative to that of stars: for power-law indiceg?| 
‘pu and ysr for BHs and stars, respectively, 2 pH & 11/4 
and 3/2 < yr < 7/4 (gi -ysr = 1) (texan ta [2017 
Long-lived stellar populations usually have A < 0.1, and the 
Galactic center is expected to be strongly segregated (Morris | 
[1993]. [Alexander & Hopman (2009) 2009) have shown the RII of 
strong mass segregation on the density distribution of a model 
stellar population around Sgr A* (mass Msmpu = 4 x 10° Mo); 
the modelled population includes main sequence stars (MS~ 1 
Mo) and stellar remnants, including white dwarfs (WD~ 0.6 
Mo), neutron stars (NS~ 1.4 Mọ) and black holes (BH~ 10 
Mo). Their study demonstrates that the heavier objects produce 
steeper density distributions via the mass segregation process. 
We note that, in the radial range 0.04-1 pc, the power-law in- 
dex of the 2D surface-density profile of the HCRs is l Hcr = 
1.6 + 0.2 (this paper) which is much greater than the corre- 
sponding value of Ix, for the K, stars in the same radial range 
(Gallego-Cano et al. 2018). The difference between the power- 
law indices (l'ucg —I'x,) is ~1. The flat power-law profile for 
K, stars provides evidence for the presence of a dynamically 


relaxed stellar cusp at the Galactic center (Gallego-Cano et al. 
2018||Schódel |2018). On the other hand, with the arguments 


presented 1n the previous paragraph, the dense group of HCRs 
reported in this paper, with its much steeper radial distribution 
(see Fig. A2) than that of the faint K, stars, could represent a 
population of massive stellar remnants that are mass-segregated 
in the nuclear star cluster at the Galactic center. The HCRs 
are associated with active massive stellar remnants having rel- 
atively higher radio luminosities. The surface density profile 
of the HCRs (see Fig. A2) may be subject to change when a 
deeper radio survey is carried out with a much lower flux den- 
sity cutoff. Nevertheless, the results from our study of HCRs in 
the central parsec provides evidence consistent with the pres- 
ence of a distribution of massive stellar remnants that is a steep 
function of the radial distance from Sgr A*. 

In addition, if the massive stellar remnants associated with 
the HCRs had migrated inward via dynamical friction, the bump 
in surface density, Xycpr, at 0.1-0.3 pc could be attributed to an 
accumulation in that radial range since the dynamical friction 
force acting on a massive object ceases at roughly half the ra- 
dius of the stellar core (~ 0.25 pc) (Merritt |2010). A maximum 
in the density profile of massive remnants 1s also predicted to 
occure at ~ 0.2 pc at a time >1 Gyr in a dynamical evolution 
model, and the bump slowly grows and migrates inward due 
to the friction produced by fast-moving stars inside these radii 


(Antonini & Merritt [2012). The Sycr-bump may serve as an 


additional observational signature of massive stellar remnants 


3 Throughout this paper, the lower-case ~y is used for power-law indices of 3D 
radial density distributions while the upper-case [ stands for power-law indices 
of 2D radial density, or the surface density, distributions. 


as a consequence of stellar dynamical processes in galactic nu- 
clei. 


5.3.2. Radiation spectrum 


Hailey et al. |(2018) recently reported the identifications of a 


dozen low-mass back hole X-ray binary candidates within the 
central parsec, implying the presence of a large population of 
X-ray binaries and isolated black holes residing within that vol- 
ume (Generozov et al. |2018). Most of them are in a quies- 
cent state. X-ray and radio flares and outbursts from these stel- 
lar remnants have been discovered over the past three decades, 


such as the magnetar SGR J1745-2900 (Kennea et al. |2013| 
2020 
12013; 2017), the microquasar of the X- 
ray transient XJ174540.0290031 
Zhao et al. 2009) and the Galactic center transient 
(GCT, see|Zhao et al. |1992) as well as the X-ray PWN candi- 
date G359.95-0.04 (Baganoff et al. [2003| 
Muno et al. [2008) which is likely powered by a neutron star. 

With a flux-density range bewteen ~0.1 and a few mJy at 33 
and 44.6 GHz, the HCRs appear to be candidate radio counter- 
parts of the old massive stellar remnants produced in the end 
of stellar evolution as expected. Most of them are in a quies- 
cent state. Although their progenitors and ages are unknown, 
the 38 steep spectrum HCRs determined from the JVLA obser- 
vation at Ka and Q bands provide important clues on the nature 
of the radio radiation, affirming nonthermal radiation with a 
steep power law in the distribution of relativistic electrons. The 
nonthermal emission could be produced by synchrotron jets or 
outflows that were launched from the compact stellar remnants 
powered by accretion from the dense, local medium. 

As described above, one of the 17 inverted-spectrum HCRs is 
the GC magnetar, SGR J1745-2900, which shows high-frequency 
pulses up to 291 GHz (Torne et al. 2017). The continuum emis- 
sion from this magnetar has been firmly detected at high radio 
frequencies up to 340 GHz and this paper). 
The inverted-spectrum of the mm-component of a magnetar to- 
wards the submillimeter appears to be a remarkable radio wave- 
length signature, as predicted by the dynamical model of 
using a persistent flow of electron-position 
plasma. The configuration of the magnetosphere of magnetars, 
created by enforced electric current and radiative drag together, 
is subject to two-stream instability. Consequently a relatively 
hard radio spectrum that is predicted to emerge, perhaps extend- 
ing to IR/optical/UV wavelengths, and is expected to be gener- 
ated because the instability leads to a large plasma density, and 
thus a large plasma frequency (Beloborodov |2013). The theory 
also predicts a large electric current associated with the radio- 
submillimeter emission from magnetars (Beloborodov |2013}, 
producing a bright radiation beam much broader than the typi- 
cal pulse width of normal pulsars with similar periods 
(2007). A valuable next step will be to use this 
theory to formulate predicted shapes of radio spectra and pulse 
profiles of magnetars for comparison with observations. 

The suggestion has been made that pulsars formed at or near 
the Galactic center might mostly to be magnetars, given the 
rather highly magnetized interstellar medium of the GC region 
that could produce highly magnetized massive stars. Subse- 
quently, strongly magnetized neutron star remnants form be- 
cause of the collapse of stellar core and concentration of the 


flux-frozen magnetic field (Dexter & O'Leary |2014 
[[2014). However, because the magnetic flux within the core 
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of a massive star as well as within a neutron star can undergo 
considerable evolution owing, for example, to dynamo action 


occurring inside the star and the neutron star, (e.g., [Duncan &| 
1992} [Thompson & Duncan |1993), this scenario 


remains rather speculative. In any case, the discovery of SGR 
J1745-2900 makes the central parsec a potentially interesting 
region to search for new magnetars. Although the recent forma- 
tion of a large number of massive stars there could have led to a 
population of neutron star remnants, if those remnants are pre- 
dominantly highly magnetized, then, as with magnetars in gen- 
eral, they would have short lifetimes as pulsars (~ 10? — 10? yrs) 
because of the powerful spindown torque associated with the 
interaction of the neutron star's magnetic field with the plasma 


in its environment (Harding et al. [1999 2011} 
Dexter & O’ Leary ;|Kaspi & Beloborodov . The 


time period during which a magnetar is an observable pulsar 
is therefore more than two orders of magnitude shorter than 
the lifetime of the massive stars observed in the young nuclear 
cluster (~2 - 10 Myrs). So we might therefore expect only a 
few (or zero) magnetars to be found as pulsars at any one time 
in the central parsec if the above speculation that massive GC 
stars produce highly magnetized remnants is correct. Indeed, 
in addition to the scatter-broadening that occurs primarily at 
longer radio wavelengths, the short lifetime of strongly mag- 
netized pulsars could be the main explanation for the rarity of 
pulsars at the Galactic center. The remaining open question is 
whether magnetars remain detectable as point radio continuum 
sources even after they have spun down to the point at which 
they can no longer be detectable as pulsars. If so, then we might 
consider that some of the HCRs are in that category. 

The inverted spectrum towards short wavelengths found in 
1E1547.0-5408 in addition to SGR J1745- 
2900 is also consistent with the prediction of spectral harden- 
ing at short radio wavelengths from the two-stream instability 
model (Beloborodov |2013). The combination of the inverted 
spectrum, high variability and high polarization — nearly 100% 
for the degree of linear polarization and 15% for circular polar- 
ization (Eatough et al. 12013) — appears to be unique to magne- 
tars. 

Of course, both the inverted and flat spectra of HCRs can 
also be interpreted as self-absorbed synchrotron emission from 
X-ray binaries in the hard state when the radiation is dominated 
by the emission from the corona of the compact object 


et al. 2011). 


6. CONCLUSION 


Following our 5.5-GHz JVLA survey of the Galactic center 
compact radio sources (GCCRs) within the radio bright zone, 
we have continued to explore Sgr A West using Ka- and Q- 
band data obtained by the JVLA in its A-configuration. At an 
angular resolution of 0.05 arcsec, we detected a dense group of 
hyper-compact («0.1") radio sources (HCRs) within the central 
parsec of the Galaxy. Based on a conservative 15o flux-density 
threshold, corresponding to 150 pJy at 33 GHz, we cataloged 
64 HCRs with their J2000 equatorial coordinates, flux densities 
at 33 and 44.6 GHz, angular sizes derived from 2D Gaussian 
fitting, and spectral index, a44 6/33. HCR49 is double. 

The surface-density distribution, Sycr(R), shows a local en- 
hancement or a density bump in the projected radial distance 
(R) range 1.5"—7" superimposed on a power-law distribution 
with an index of I' = 1.6+0.2. The steeper profile of the HCRs 
relative to that of the nuclear stellar cluster might result from 
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the concentration of massive stellar remnants by mass segrega- 
tion. 

The 65 HCRs divide into three spectral sub-types: 38 steep- 
spectrum (0446/33 < —0.5), 10 flat-spectrum (—0.5 < 044,6/33 < 
0.2), and 17 inverted-spectrum sources (044,6/53 > 0.2). Our sta- 
tistical analysis shows that the distribution of the steep-spectrum 
HCRs in Log(S[mJy] can be fitted to a Gaussian with a mean 
of  =—0.35 (corresponding to 0.45 mJy) and a standard devi- 
ation of ø = 0.22. We suggest that the steep-spectrum HCRs 
be regarded as candidates for a population of stellar remnants 
acting as nonthermal compact radio sources powered by accre- 
tion onto neutron stars and stellar-mass black holes, with the 
accreted matter supplied either by a binary companion or by a 
dense portion of the interstellar medium. 

The inverted-spectrum HCRs show a rising spectrum towards 
high frequencies. Five of the 17 inverted-spectrum HCRs have 
compact IR counterparts, suggesting that they are associated 
with the ionized stellar wind outflows from hot, massive stars. 
A portion of the inverted-spectrum HCRs may consist of X-ray 
binaries in the hard state, when the self-absorbed synchrotron 
emission is dominated by the corona of the compact object. 
The GC magnetar, SGR J1745-2900, belongs to the inverted- 
spectrum sub-type. Based on our analysis of 91 flux-density 
measurements of SGR J1745-2900 observed with the JVLA, 
ALMA and ATCA, we find that two distinguishable spectral 
components contribute to the averaged spectrum, separated at 
the transition frequency v, ~ 30 GHz. The cm-component is fit- 
ted to a power-law with a steep spectral index ag, =—1.5+0.6 
while the mm-component shows the inverted spectrum o4 = 
1.1 3:0.2. Our consolidation of the spectrum from the inter- 
ferometer array data is in good agreement with earlier results 
based on single-dish observations (Torne et al. [2017]. 

In addition, we reduced 225-GHz ALMA data observed with 
a FWHM beam of 0.024" x 0.017" on ten individual days within 
two weeks in October, 2017. SGR J1745-2900 was detected 
with a signal-to-noise ratio of ~70 at flux densities varying 
widely between 0.79 and 5.60 mJy on day-to-day timescales. 
An index of relative variability ZV ~ 6 is found for SGR 
J1745-2900, which compares with ZY zz 0.2 for the more 
moderately variable Sgr A*, derived from the same ALMA ob- 
servations. 

Collecting a total of 161 individual flux-density measure- 
ments of SGR J1745-2900 at radio wavelengths from prior pub- 
lished literature and this paper, we binned data into 100-day 
bins along the time axis for both cm- and mm-components sepa- 


rately, starting from the outburst on April 23, 2013 (MJD 56406). 


The radio light curves with bin-averaged flux density are com- 
pared with the X-ray light curve observed with Chandra 
[al. [2020). Except for the appreciable short-timescale variability 
at millimeter wavelengths, the long-term trend at both centime- 
ter and millimeter wavelengths is a slow, but erratic, decrease 
in radio power, in contrast to the smoothly decreasing trend in 
X-ray luminosity. 

Because many HCRs are candidates for being remnants of 
massive stars, we considered the origin and fate of such rem- 
nants, and speculated on the possible reasons for the difficulty 
of finding pulsars there. Neutron stars formed within the rela- 
tively highly magnetized central molecular zone of the Galaxy 
could themselves inherit sufficient magnetic flux that they are 
mostly born as magnetars. In that case, their ultra-strong mag- 
netic fields (10!^ - 10/6 G), acting on plasma trapped in their 
magnetospheres, apply a spin-down torque that causes them to 


have relatively short lifetimes as detectable pulsars (10? — 10? 
yr) compared to normal pulsars. 

If the inverted spectrum towards high frequencies found in 
SGR J1745-2900 and 1E1547.0-5408 is a common characteris- 
tic of active magnetars, the high-frequency bands (Ka and Q) of 
the JVLA and all the ALMA and SMA bands open a practical 
window for studying magnetars. Also, dynamical information 
on the HCRs can be acquired with high-precision proper motion 
data that can be obtained with repeated JVLA and ALMA ob- 
servations. Such measurements will supply strong constraints 
on the candidacy of the HCRs as massive stellar remnants at the 
Galactic center. 
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APPENDIX A 
DENSITY DISTRIBUTION OF HCRS AND GCCRS 


The stellar density distribution in the Galactic center can be 


fitted by two power laws (e.g., 2003 
2016), with a flat inner component and a steeper component at 


large projected radii (Schódel 12014/2018). Usually, the pro- 


jected surface density can be described by the “Nuker’ profile 


with two slopes (Lauer et al. |1995): 


RAF R \ TOA 
X(R) = X(Ry26-/o (2) L (s) | (Al) 
R Rp 
where X(R,) is the surface density at the break radius R, that 
divides the profile into two power laws. At a small radius R < 
Ry, X(R) ~ RT, usually describing the cusp in the distribution 
of old stellar population. For R >> Rp, X(R) ~ R®, fitting to the 
outer power law, which is usually steeper than that of the cusp 
(Lauer et al. 1995). The break radius R, corresponds both to 
the point at which the slope is the mean of ( and T and to the 
radius of maximum curvature of the distribution in the log;9(R)- 
logio((Rp)) coordinate system. 

In this section, we present the algorithm that we developed 
to derive the distribution of the surface density (34,) for the 
compact radio sources (CRS) as function of the projected radial 
distance from Sgr A*. 


A.1. Algorithm 


First, the radio bright zone (RBZ) is divided into 20 annular 
zones (AZ), each of them bounded by two rings at outer and 
inner radii, R(i) and R(i+ 1): 


RG) =2°"'Repz, (A2) 


where i=0,1,2...,19 and 20 and Rggz is the outer radius of the 
radio bright zone (RBZ). Fig. Al shows the 20 annular zones 
along with the distribution of GCCRs and HCRs in the RBZ 
(Fig. Ala), Sgr A West (Fig. Alb) and the central parsec (Fig. 
Alc). The radius of the outer ring, R(0) = Rrpz, corresponds 
to the largest radius within which we searched the GCCRs with 
a 10-c cutoff due to the JVLA primary beam at 5.5 GHz. The 
radius of the most inner ring, R(20) = 0.55” ~ Opwun, is ap- 
proximately the size of FWHM beam at 5.5 GHz. The zone 
AZ-10 marks the boundary, within which we found 64 HCRs 
with the JVLA A-array observations at 33 and 44.6 GHz (this 
paper). Outside this boundary, 110 GCCRs are found from the 


JVLA A-array observations at 5.5 GHz (Zhao, Morris & Goss 
2020). 


Then, we counted the CRS in each of the AZs by creating a 
grid function crs(j, k) for each of the CRSs, where k is the an- 
nular zone ID from 1 to 20 and j is the CRS source ID within an 
AZ. For given a specific grid cell (j,k), a unit value is assigned, 
i.e. crs(j, k) = 1. Given a k, the AZ-ID, the number of sources 
in an annular zone is then derived by summing over j: 


Nk) 2 V crs(j,k). (A3) 
j 


In addition, we also created grid functions of s( j,k) and o(j,k) 
for flux-density and uncertainty, respectively. Corresponding 
to a specific source ID: (j,k), the values of flux density and 
l-o uncertainty of the CRS are assigned to s(j,k) and o(j,k). 
The accumulated flux density of CRS in an annular zone can be 
determined: 


N(k) 


Sik) = Y ^s), (A4) 


i=l 


and 


N(k) 
o(k)= |) Lo ky. (A5) 
j=l 


We note that, for the GCCRs, s(j,k) and o(j,k) correspond 
to the mean and 1-o error of the mean, which are determined 
with the measurements from the three-epoch observations at 5.5 


GHz (Zhao, Morris & Goss |2020). For the HCRs, we used Ka- 
band flux density and 1-c error at 33 GHz (this paper). The 


area for the k-th AZ can be computed: 


A(k) = x [R - R(c- 17] . (A6) 
Furthermore, we also estimated a flux-density-based CRS counts 
Ns(k): 
S(k) 


Ns(k) = 
st ) S,,5.5 GHz’ 


(A7) 


fork < 11 and 


Nso = 0. (A8) 
1,33 GHz 
for k > 11. Here $,,55 GHz and S,,33 GHz are the mean flux- 
density values for the GCCRs at 5.5 GHz and HCRs at 33 GHz. 
Given the sensitivity cutoff and the variability of the CRSs, 
the direct source counts N(k) represents a lower limit of the 
k-th AZ. The flux-density-based CRS counts Ns(k) is usually 
greater than the direct source counts, i.e., Ns(k) > N(k). Thus, 
using the averaged values of [Ns(k)+N(k)] /2 we derived the 
surface density of CRSs as a function of radius from Sgr A*: 
oo (A9) 
2A(k) 
The uncertainty of surface density is estimated as well: 


[Ns(k) -N(K)] 
As,,(k) = 2A0) 
where k = 1,2,3, ... and 20. Thus, the upper and lower limits 
in surface density correspond to the flux-density-based density 
(Ns(k)/A(K)) and the actual-source-counted density (N(K) / A(K)), 
respectively. 
We then calculated the mean radius given the k-th annular 
zone: 


(A10) 


R(k)+R(K+1 
Raz(k) = = 


and the uncertainty, or more properly half of the annular zone 
width: 


(A11) 


Ara (k) = ie AT D ; (A12) 
where R(k) is defined in Eq.(A1) and k = 1,2,3, ... and 20. 

Fig. A2 shows the distribution of the surface density of CRSs 
Mas (R) as a function of radial distance (R) from Sgr A*. The 
black symbols represent the data at a large projected radius 
(R > 13"), determined from the JVLA observations at 5.5 GHz 
(Zhao, Morris & Goss [2020) while the symbols for the inner 
annular zones (R < 13") mark the data derived from the Ka- 
band observations at 33 GHz. We note that the CRS counts, as 
well as the surface density 3, (light-grey symbols in Fig. A2), 
may be underestimated within the inner annular zone (AZ=11 - 
20) due to the fact that the data HCRs observed at 33 GHz cor- 
responds to a higher effective cutoff limit at 5.5 GHz since most 
of the HCRs have a steep spectrum. With an averaged spectral 
index aycrs = —0.42 + 0.12 using all 64 HCRs' data, we con- 
verted the 33-GHz flux densities to 5.5-GHz. The equivalent 
cutoff limit for the inner AZs is about twice as high as the limit 
searched for GCCRs in the RBZ (Zhao, Morris & Goss [2020). 
So, we multiplied the surface density of the inner AZ(k>11) by 
a correction factor: 


n = (5.5/33) CR (=) x24, (A13) 
05.5 

assuming that the HCR counts are inversely proportional to the 

cutoff value in flux density. The corrected surface density is 

then computed: 


Ms Dos (A14) 
as indicated with the red symbols which appear to be aligned 
better with the surface density data of the outer radii obtained 
from the observations at 5.5 GHz. 
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Fic. Al.— The radio bright zone (RBZ) is divided into 20 annular zones 
(AZ) along the radial direction from Sgr A*. Spatial distribution of GC com- 
pact radio sources or GCCRs (light blue symbols, crosses for variables and 
dots for non-variables) in annular zones is displayed over the RBZ (AZ 1 - 6) 
(top panel: Fig. Ala) and SgrA West (AZ 7-11) (middle panel: Fig. Alb). 
The HCRs (light blue dots) are distributed in the central parsec (AZ 12-20) 
(bottom panel: Fig. Alc). The thick red ring marks the 11-th ring (Eq. A2), 
R(11) = 13", ~ 0.5 pc. In Figs. Ala and Alb, the symbol of a black cross-dot 


marks the Galactic center transient, the GCT (Zhao et al. 1992). The magnetar 
Ac. 


J1745-2900 and IRS 13E are labelled in Fig. 


A.2. Power-law fit to the Xor distribution 


The surface-density 3X, data at large radii (from AZO1 to 
AZ09) can be described by a power law while the data at small 
radii (from AZ10 to AZ18) show a bump (%,,,-bump), indicat- 
ing a locally enhanced source density 3, at the radial distance 
range of R = 1.5"—7.0" (from AZ-17 to AZ-14) superimposed 
on a slope. We made a linear regression fit to a function: 


Log o(%ers) = agccr — Logio(&), (A15) 


using the outer nine AZ's data points (GCCRs) with a weight 
of wt(i) — AS. (k). The best fit derived from the least squares 
regression (LSR), as indicated by a straight blue line in Fig. 
A2, gives a power law function X, cx R^? (R > 13") with B= 
2.0+0.1 and AGCCR = 4.3+0.2. 

For the inner eight annular-ring data (HCRs), we excluded 
the four data points that are involved in the Xers-bump (from 
AZ14 to AZI7) and fit the remaining four data points with a 
linear function: 


Logio(3as) = ancr — VLogio(R). (A16) 


Based on a least-squares fitting, we derived X, cx RT (R< 
13") with I = 1.6 + 0.2 and aycr = 3.7 4: 0.2, as indicated with 
the straight orange line in the logarithmic plot (Fig. A2a). The 
Xcrs-bump is characterized by the high surface-density zones 
AZI15 and AZ16; the former contains the GC magnetar J1745- 
2900 and the latter is associated with IRS 13E. We note that the 
surface brightness in the diffuse Bry narrow band filter image 
(Fig. 6a of |Schódel |(2018)) also shows a bump in the range of 
radius between R = 1.5" and 7.0", suggesting that the counts of 
CRSs in the radial range may be partially caused by the local 
hyper compact HII sources associated with younger massive 
stellar objects. However, excluding the flat-spectrum HCRs, the 
Xer -bump is also present (see Fig. A2b). We further checked 
the distribution only with the steep-spectrum HCRs and found 
that the 34,,—bump is still present in the radial range between 
Rz 1.5" and 7.0" (see Fig. A2c). 

Due to a small number of data points and at least four free 
parameters in the “Nuker” model of Eq(A1), we made no at- 
tempt to fit the entire curve for deriving the break radius Rp. 
However, from the intersection between the two straight lines 
of Eq(A15) and Eq(A16), we found an approximate break ra- 
dius R, ~ 30", or 1.1 pc, which appears to correspond to the 
radius of the ionized region. Therefore, in the main text, we use 
R, = 30" as the break radius that separates the two distinctive 
power-law distributions. 
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